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sThi*  report  describes  work  performed M.I.T.  Lincoln  Laboratory,  between  I  October  1979  and  I  March 
1983,  C®  evaluate  the  use  of  the  Global  Positioning  System  (CPS)  for  low-cost  civil  air  navigation. 
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The  report  describes  a  GPS  Test  and  EvaluatlorJjSyslem  developed  jointly  by  M.I.T.  Lincoln  -Laboratory^’  / 
Stanford  Telecommunications,  Inc.,  and  Intermetrics,  Inc.,  using  techniques  that  could  lead  to  low-cost  com¬ 
mercial  avionics.  System  performance  results  obtained  in  the  laboratory  and  during  flight  tests  are  provided 
which  demonstrate  compliance  with  current  and  future  navigation  accuracy  requirements  for  enroute,  termi¬ 
nal  and  non-precision  approach  flight  paths.  The  report  also  includes  functional  specifications  for  a  low-cost 
GPS  navigation  system  for  civil  aircraft.  Tkt  C,  P  5  7V  £~  — -y 

The  GPS  Test  and  Evaluation  (ystemdesign  was  based  on  two  important  features:^)  automatic  tracking  of 
all  visible  satellites  (rather  than  a  minimum  set  of  four)jand(2)  a  dual-channel  GPS  C/A  code  receiver.  Track¬ 
ing  all  visible  satellites  allows  the  system  to  maintain  continuous  navigation  when  a  satellite  sets  or  is  momen¬ 
tarily  masked  during  aircraft  maneuvers.  The  dual-channel  receiver  dedicates  one  channel  to  pseudo-range 
measurements,  and  the  other  channel  to  acquiring  new  satellites  as  they  become  visible.  These  two  features, 
validated  by  flight  test,  allow  the  system  to  provide  continuous  navigation  updates  during  critical  aircraft 
maneuvers,  such  as  non-precision  approaches,  and  during  satellite  constellation  changes.  Of|0/  vo  A-  r  -  S  u 
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1.0  INTRODUCTION  AND  SUMMARY 


This  report  describes  the  work  performed  by  M.I.T.  Lincoln  Laboratory, 
between  1  October  1979  and  1  March  1983  to  evaluate  the  use  of  the  Global 
Positioning  System  (GPS)  for  low-cost  civil  air  navigation.  The  effort  was 
supported  by  the  Federal  Aviation  Administration  through  Interagency  Agreement 
D0T-FA79WAI-091  between  the  FAA  and  the  United  States  Air  Force. 

M.I.T.  Lincoln  Laboratory  was  tasked  by  the  FAA  to  develop  a  GPS  Test  and 
Evaluation  System  and  to  flight  test  the  system  in  a  typical  general  aviation 
aircraft.  The  system  was  required  to  meet  FAA  accuracy  requirements  for  area 
navigation  systems  and  to  use  techniques  that  could  lead  to  low-cost 
commercial  avionics.  This  report  describes  the  system  design,  and  provides 
system  performance  results  from  laboratory  and  flight  tests.  The  report  also 
includes  functional  specifications  for  a  low-cost  GPS  navigation  system  based 
on  validated  design  concepts. 

The  GPS  Test  and  Evaluation  system  design  was  based  on  two  important 
features:  1)  automatic  tracking  of  all  visible  satellites  (rather  than  a 
minimum  set  of  four)  and  2)  a  dual-channel  GPS  C/A  code  receiver.  Tracking 
all  visible  satellites  allows  the  system  to  maintain  continuous  navigation 
when  a  satellite  sets  or  is  momentarily  masked  during  aircraft  maneuvers.  The 
dual-channel  receiver  dedicates  one  channel  to  pseudo-range  measurements,  and 
the  other  channel  to  acquiring  new  satellites  as  they  become  visible.  These 
two  features  allow  the  system  to  provide  continuous  navigation  updates  during 
aircraft  maneuvers  such  as  non-precision  approaches  and  during  satellite 
constellation  changes. 

The  major  elements  of  the  GPS  Test  and  Evaluation  System  are: 

•  Dual  Channel  Receiver  -  a  unit  developed  by  Stanford 
Telecommunications,  Inc.  that  sequentially  tracks  satellites  on  one 
channel  while  the  other  channel  acquires  new  satellites  and 
demodulates  data;  provides  pseudo-range  estimates  and  satellite  data 
to  the  Position  Processor. 

•  Position  Processor  -  a  microcomputer  with  software  developed  by 
Intermetrics ,  Inc.  that  manages  the  receiver  channels  and  computes 
geodetic  position  estimates. 

•  Navigation  Processor  -  a  microcomputer  that  translates  geodetic 
estimates  into  displays  suitable  for  pilot  navigation. 

•  Pilot  Displays  -  a  standard  Course  Deviation  Indicator/Omni -Bearing 
Selector  (CDI/OBS)  and  a  Control  and  Display  Unit  (CDU),  designed  to 
provide  navigation  data  to  the  pilot  in  a  format  consistent  with 
current  civil  navigation  practices. 

The  system  also  includes  extensive  instrumentation  and  data  recording 
capabilities  for  post-flight  data  analysis. 


The  performance  of  the  GPS  Test  and  Evaluation  System  was  measured  during 
laboratory  tests  and  during  flight  tests  in  a  typical  twin-engine  general 
aviation  aircraft,  a  Rockwell  Aero  Commander.  The  laboratory  test  showed  that 
the  typical  static  RMS  horizontal  position  error  was  93  feet.  The  flight 
tests  showed  a  RMS  horizontal  error  of  189  feet  and  a  95%  confidence 
horizontal  error  of  333  feet  in  typical  general  aviation  operations.  As  shown 
below,  the  measured  accuracy  was  well  within  the  requirements  of  FAA  Advisory 
Circular  90-45A  for  two-dimensional  area  navigation  and  essentially  met  the 
proposed  328  ft.,  95%  Federal  Navigation  Plan  requirement  for  non-precision 
approach. 

Navigation  Accuracy  Summary 

Horizontal  (95%) 

un 

FAA  Advisory  Circular  90-45A 

Approach,  ascent/descent  1824. 

Federal  Radionavigation  Plan 

Non-precision  approach  328. 

Accuracy  Achieved 

Level  and  turning  mix  333. 

Operation  of  the  GPS  Test  and  Evaluation  System  was  evaluated  in  areas  of 
mountainous  terrain,  at  a  large  urban  airport  and  in  typical  general  aviation 
operations.  The  system  was  shorn  to  provide  continuous  navigation  service 
during  30°  bank-angle  turns  and  was  able  to  track  satellites  with  elevation 
angles  as  low  as  5°.  Finally,  the  system  appeared  to  be  compatible  with 
existing  air-traffic  control  procedures  and  air-crew  practices. 

Although  the  experimental  GPS  was  found  to  perform  well  as  a  practical 
general  aviation  navigation  system,  it  was  physically  large.  This  was 
necessary  in  order  that  the  system  be  constructed  from  readily  available 
components  and  that  it  provide  for  experimental  flexibility.  However,  the 
increasing  use  of  VLSI  digital  techniques  is  likely  to  result  in  a  system 
design  which  can  be  produced  in  an  avionics  package  comparable  in  size  to  that 
of  present  day  commercially  available  area  navigation  equipment.  The  cost,  as 
estimated  by  ARINC  Research,  Inc.,  of  a  GPS  navigator  based  on  the  Lincoln 
design  was  $8500  in  1982  dollars  assuming  the  use  of  circa-1990 
integrated-circuit  technology. 


The  remainder  of  this  report  is  divided  into  four  sections.  Section  2.0 
provides  an  overview  of  the  project  and  details  the  performance  requirements 
and  goals.  Section  3.0  describes  the  system  design,  including  the  hardware 
characteristics  and  software  architecture.  Section  4.0  provides  the  system 
performance  results,  including  ground  and  flight  test  measurements. 

Section  5.0  gives  functional  requirements  of  a  general  aviation  GPS  navigation 
system  as  determined  from  the  results  of  this  project. 


2.0  PROJECT  OVERVIEW 


2. 1  Background  and  Objectives 

The  NAVSTAR  Global  Positioning  System  (GPS)  is  a  satellite-based  system 
currently  being  developed  by  the  Department  of  Defense  to  provide  position 
information  to  suitably  equipped  military  users.  Such  users  will  receive 
world-wide,  continuous,  real-time,  all  weather,  precision  navigation  data  from 
a  constellation  of  18  satellites  in  12  hour  orbits  (Ref.  1).  Six  satellites 
have  been  placed  in  orbit  and  ground  support  facilities  have  been  developed  to 
support  user  equipment  development. 

Previous  studies  (Refs.  2,  3)  suggested  that  a  GPS  navigator  meeting  FAA 
performance  requirements  could  be  developed  at  a  cost  consistent  with  general 
aviation  use.  However,  the  performance  of  low  cost  GPS  receiver  designs  had 
not  been  fully  determined  by  field  measurement.  To  fill  this  gap,  M.I.T. 
Lincoln  Laboratory  was  tasked  by  the  Federal  Aviation  Adminstration  to 

a)  develop  GPS  Test  and  Evaluation  Equipment  functionally  consistent  with  a 
low  cost  design  and  FAA  two-dimensional  (2-D)  navigation  requirements, 

b)  evaluate  the  design  by  field  measurement,  and  c)  develop  a  functional 
specification  for  a  general  aviation  GPS  receiver.  Lincoln  Laboratory 
contracted  with  Stanford  Telecommunications,  Inc.,  to  develop  a  dual  channel 
receiver,  and  Intermetrlcs ,  Inc.,  to  develop  receiver  management  and  position 
estimation  software.  To  establish  cost,  the  FAA  tasked  ARINC  Research,  Inc., 
to  estimate  the  cost  of  a  commercial  version  of  the  design  developed  at 
Lincoln  Laboratory;  the  results  of  that  cost  study  have  been  reported  by  ARINC 
Research  in  a  separate  document  (Ref.  4). 

2.2  Current  Civil  Navigation  Requirements 

In  order  that  a  navigation  system  be  acceptable  for  civil  aviation  it 
must  satisfy  needs  within  four  major  areas  of  concern:  accuracy,  reliability 
integrity,  and  compatibility. 

2.2.1  Accuracy 

The  basic  output  of  the  typical  GPS  navigation  system  is  a  geodetic 
(latitude-longitude)  position  estimate.  This  is  to  be  contrasted  with  low 
cost  VOR/DME  navigation  receivers  which  produce  rho-theta  (range-bearing) 
estimates  with  respect  to  a  ground  facility.  The  relevant  specifications  from 
FAA  Advisory  Circular  90-45A,  (Ref.  5)  for  GPS  are  shown  summarized  in 
Table  2-1. 

Position  accuracy  should  not  be  considered  apart  from  update  rate. 
Equipment  intended  for  use  in  non-precision  approaches  should  provide 
navigation  information  that  is  essentially  continuous  with  interruptions  no 
longer  than  would  result  from  switching  from  one  pre-programmed  waypoint  to 
another.  Also,  navigation  service  should  be  continuously  available  during 
flight  in  any  directlon/cllmb/descend  profile  approved  by  the  aircraft 
manufacturer.  The  navigation  service  shall,  per  AC90-45A,  be  restored  (if 
temporarily  lost  due  to  bank-induced  fades)  within  5  seconds  after  the 
completion  of  any  allowable  maneuver.  Also,  the  time  lag  between  data 
measurements  and  displayed  position  must  not  be  operationally  significant. 


TABLE  2-1 


AC90-45A  2-D  NAVIGATION  REQUIREMENTS 


Position  Accuracy  Requirements  -  95%  Conf. 
(Averaged  Over  One  Update  Cycle) 

Cross-Track 

Along-Track 

Flight  Technical  Error1 2  (FTE) 

Enroute^ 

2.0  nml 

— 

Terminal3 * 5 

1.0  nml 

— 

Non-Preclslon  Approach* 

0.5  nml 

— 

Navigation  Equipment  Error 

Enroute 

1.5  nml 

1.5  nml 

Terminal 

1.1  nml 

1.1  nml 

Non-Preclslon  Approach 

0.3  nml 

0.3  nml 

Total  Error3 

Enroute 

2.5  nml 

1.5  nml 

Terminal 

1.5  nml 

1.1  nml 

Non-Preclslon  Approach 

0.6  nml 

0.3  nml 

1.  FTE  accounts  for  deviations  due  to  display  Interpretation,  pilot 
response  and  aircraft  response. 

2.  Knroute  «  Cruising  flight  between  terminal  areas. 

3.  Terminal  »  Between  enroute  and  approach;  normally  below  18000  feet  and 

within  SO  miles  of  the  airport. 

A.  Non-Preclslon  Approach  ■  Between  final  approach  waypoint  and  airport. 

5.  Total  Error  ■  ESS  combination  of  FTE  and  navigation  equipment  errors. 
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A a  part  of  the  test  and  evaluation  equipment  design  development ,  the 
accuracy  and  update  specifications  in  Table  2-2  were  established.  The  dynamic 
requirement  is  the  Non-Precision  Approach  approach  tolerance  of  AC90-45A.  The 
static  requirement  is  derived  from  a  generic  low-cost  receiver  model.  An 
additional  requirement,  Time-to-First-Fix,  was  established  based  on  the 
typical  time  between  first  turn-on  of  the  avionics  and  when  the  general 
aviation  aircraft  is  positioned  near  the  departure  runway. 

2.2.2  Reliability 

The  CPS  navigation  system  encompassing  the  user  equipment,  satellite 
vehicles  and  ground  support  equipment,  must  have  a  combined  reliability  equal 
to  or  surpassing  that  of  alternative  navigation  systems.  The  GPS  user 
equipment  must  therefore  provide  navigation  data  without  operationally 
significant  outages  due  to  fades  or  multipath,  assuming  that  the  NAVSTAR 
constellation  provides  acceptable  coverage. 

2.2.3  Integrity 

The  GPS  navigation  system  must  not  provide  misleading  information  under 
any  conditions  of  operational  significance.  It  is  therefore  necessary  that 
the  GPS  receiver  continually  monitor  its  own  performance  and  indicate  to  the 
pilot  when  the  navigation  Information  displayed  is  no  longer  in  compliance 
with  the  accuracy  requirements. 

Further,  provisions  must  be  Incorporated  to  allow  the  pilot  to  verify 
that  the  navigation  information  is  accurate  using  either  built-in  test 
equipment,  an  auxiliary  test  system,  or  a  procedural  check. 

2.2.4  Compatibility 

The  GPS  navigation  system  must  provide  a  pilot  interface  which  is 
compatible  with  existing  air  navigation  systems.  This  requirement  stems  from 
the  way  that  pilots  are  accustomed  to  using  aircraft  navigation  systems. 


2.3  Future  Civil  Navigation  Requirements 

To  anticipate  the  more  stringent  navigation  requirements  likely  in  the 
future,  the  accuracy  and  update  rate  goals  in  Table  2-3  were  established  for 
the  GPS  Test  and  Evaluation  (T  and  B)  equipment. 

Later,  after  the  GPS  T  and  E  Equipment  specifications  were  frosen  in 
procurement  specifications,  the  Federal  Radio  Navigation  Plan  (FRP)  (Ref.  6) 
was  completed  and  published.  The  FRP  established  new  future  navigation 
requirements  that  are  summarised  in  Table  2-4. 

The  rationale  for  increasing  tha  accuracy  is  to  provide  a  service 
equivalent  to  that  provided  by  on-airport  VOR's  which  now  exist  at 
approximately  302  of  all  U.S.  airfields. 


TABLE  2-2 


GPS  TEST  AND  EVALUATION  EQUIPMENT  ACCURACY  AND  UPDATE  REQUIREMENTS 


Horizontal  Dynamic  Position  Error* 

0.3  NM,  95Z 

Horizontal  Static  Position  Error^ 

335  ft.,  95Z 

Update  Rate 

0.5  Hz,  min 

Time  to  First  Fix3 

Assumptions: 

Velocity  <  200  KTS 

Acceleration  <  0.5g 

LI  C/A  only 

6  minutes,  max 

No  artlflcal  degradation  of  the  GPS  navigation  message  data 
or  Intentional  electromagnetic  Interference 

1.  For  HDOP*  <  5,  bank  angle  >  30°. 

2.  For  HDOP  <2.0  and  a  2  minute  Integration. 

3.  From  power-on  until  cockpit  display  of  position  fix. 

horizontal  Dilution  of  Precision. 
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Horizontal  Position  Error* 
Level  Flight 
Turning  Flight 


500  ft.,  95% 
1000  ft.,  95% 


Vertical  Position  Error^ 
Enroute  level  flight 
Terminal  level  flight 
Approach  level  flight 


163  ft.,  95% 
163  ft.,  95% 
88  ft.,  95% 


Enroute  Ascend /Descent 
Terminal  Ascend /Descent 
Approach  Ascend /Descent 


275  ft.,  95% 
205  ft.,  95% 
116  ft.,  95% 


Update  Rate 


1  Hz,  min 


1.  For  HD0P  <  5,  bank  angle  <  30°,  ground  speed  <  200  KTS. 


2.  From  AC90-45A  (Total  error  less  FTE,  Table  B,  Appendix  A) 
conditions  as  Horizontal  Goals. 


Phase 


Altitude 


Horizontal  Position  Error 


Enroute 

500  ft.  to  FL  180* 

1000  meters 
(3280  ft.) 

(2  drms) 

Terminal 

500  ft.  to  FL  180 

500  meters 
(1650  ft.) 

(95Z) 

Non-Preelslon 

Approach 

*FL  180  means 

250  to  3000  ft.  above  surface 

flight  level  18,000  feet. 

100  meters 
(328  ft.) 

(95Z) 

for  the  development  and  testing  of 


2.4  Aasumgtlona 

The  following  assumptions  were  made 
the  GPS  T  and  S  Equipment : 

a.  The  satellites  will  radiata  the  LI  signal  such  that  the  received 
signal  at  the  civil  users  GPS  antenna  Input  will  aeet  or  exceed  the 
levels  specified  In  the  space  segnent  User  Interface  ICD,  (Ref.  7). 

b.  The  radiated  signals  from  the  GPS  satellites  will  not  be  artlfically 
degraded  to  reduce  accuracy. 

c.  The  T  and  E  Equipment  will  not  be  exposed  to  Intentional  radio 
frequency  Interference. 

2.5  General  Approach 
2*5.1  Development 

The  development  of  the  GPS  T  and  E  Equipment  was.  In  part,  based  on  past 
work  (Ref.  7)  and  on  the  DOD  GPS  user  equipment  developments  that  suggested: 

a.  Designs  which  tracked  all  satellites  In  view  (a  significant 
advantage)  could  be  realised  at  reasonable  cost. 

b.  Designs  which  convert  functions  from  hardware  to  firmware  using 
microcomputers  would  reduce  cost,  and 

c.  An  LI  C/A-only  receiver  would  be  able  to  provide  acceptable  accuracy 
for  civil  general  aviation  users. 

The  need  for  continuous  navigation  service,  especially  during 
non-precision  approaches,  led  to  the  requirement  that  the  T  and  E  Equipment 
contain  two  independent  channels.  The  need  to  evaluate  performance  under  a 
variety  of  stressing  situations  resulted  in  the  requirement  for  extensive 
Instrumentation.  The  requirement  to  evaluate  operational  Issues  lead  to  the 
definition  of  a  general  purpose  cockpit  display  consistent  with  current  civil 
navigation  procedures. 

It  was  assumed  that  during  the  Initial  general  aviation  use  of  GPS,  the 
pilot  Interface  would  be  consistent  with  the  VOR/DKE  system,  with  RHAV 
features  as  added  options.  This  Implied  that  the  following  features  would  be 
available  on  Initial  GPS  receivers  and  therefore  should  be  part  of  the  T  and  E 
Equipment: 

a.  Automatic  Operation.  The  receiver  should  automatically  acquire  and 
track  satellites,  estimate  position  and  continually  monitor 
performance  without  operator  Intervention. 

b.  Course  Deviation  Indicator  (CPI)  Interface.  The  receiver  should 
accept  desired  course  Waring  From  an  Omni  Bearing  Selector  (OBS)  and 
provide  signals  to  drive  the  Course-Deviation  Indicator  (CDI). 


c. 


Magnetic  B— ring.  A  np  should  bo  provided  to  apply  local  nagnetlc 
variation  corrections  (to  within  1  1  degree);  the  nap  to  be  stored  in 
a  non-volatile  memory. 

d.  Way  Points.  The  systea  should  accept  way  points  in  the  fora  of 
unaablguous  designators  such  as  the  current  letter  codes.  The 
receiver  could  then  aalntaln  several  thousand  waypoints  with 
their  associated  geographical  coordinates  in  a  non-volatile  aaaory. 

The  concern  for  integrity  led  to  the  inclusion  of  performance  monitoring 
and  fall-soft  techniques.  Testa  such  as  satellite  health  checks,  pseudo-range 
consistency  tests,  position  estimate  variance  tests  and  geometry  tests  were 
provided. 

Finally,  low  cost  was  emphasised  by  designing  to  requireasnts  that  are 
consistent  with  civil  aviation  navigation  and  by  using  technology  which  can  be 
expected  to  become  inexpensive  in  the  current  decade. 

The  design  and  construction  of  the  equipaent  was  accomplished  during 
1980-82.  As  mentioned  previously.  Standard  Telecoamunicatlona,  Inc., 
developed  the  dual  channel  receiver  and  Interaetrlca,  Inc.,  developed  software 
which  managed  the  receivers  and  eatlaated  poaltlon.  Lincoln  Laboratory 
developed  the  navigation  aoftware,  lnatrumentatlon.  Control  and  Display  Doit 
(CDO),  and  ground  support  facilities,  conducted  ground  and  flight 
aeaaureaanta ,  and  evaluated  systea  performance. 

2.5.2  Evaluation 


In  keeping  with  the  principal  concerns,  field  evaluation  of  the 
GPS  T  and  E  equipment  focused  on  a)  link  margins  and  position  estimation 
techniques,  b)  performance  monitoring  and  fall-soft  techniques,  and 

c)  operational  performance  in  typical  general  aviation  flight  operations  with 
emphasis  on  non-precision  approaches. 

Link  margins  were  evaluated  by  determining  the  effecta  of: 

a)  the  number  and  location  of  available  satellites, 

b)  received  signal  strength, 

c)  multipath  and  terrain  blockage, 

d)  electromagnetic  interference  (EMI), 

e)  antenna  shielding,  and 

f)  receiver  acquisition  and  tracking 

on  the  acqulaltlon  time,  positional  accuracy  and  update  reliability.  Accuracy 
was  measured  using  truth  derived  from  a  ground  tracker  while  flying  the 
T  and  E  Equipment  in  a  general  aviation  aircraft.  The  position  estimation 
technique,  including  ita  smoothing  filter,  was  evaluated  to  assess  the  effect 
of  signal  margin  and  aircraft  dynamics  on  the  position  estimate  variance. 


wtlBitd  within  the 
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3.0  OPS  TEST  AMD  EVALUATION  EQUIPMENT  DESIGN 


3.1  Syitei  Design 

The  OPS  T  and  E  Equipment  is  shown  in  Fig.  3-1.  This  design  provides 
Ugh  rate  (5  He)  sequential  processing  of  all  visible  GPS  satellites  in  order 
to  support  a  position  estiaate  update  rate  of  1.0  Hs.  Two  channels  are 
provided  in  order  for  one  channel  to  be  dedicated  to  Ugh  rate  sequential 
pseudo-range  processing  while  the  other  is  available  to  deaodulate  satellite 
navigation  messages.  This  assures  continuous  navigation  service  during 
non-precision  approaches. 

The  key  features  of  the  systea  design  are: 

a.  Autoaatlc  Operation.  The  receiver  autoaatically  acquires  all 
satellites  in  view  and  provides  a  first  fix  to  the  pilot  in  less  than 
6  minutes.  If  required  (after  being  powered-off  for  aore  than  30 
days),  it  autoaatically  acquires  fresh  alaanac  data,  which  Increases 
the  tiae-to-flrst  fix  (TTPF)  to  16  nloutes.  The  receiver  noraally 
requires  no  operator  Intervention  during  use  other  than  such 
navigation  functions  as  waypoint  entry. 

b.  Perforaance  Monitoring .  The  receiver  continually  aonitors  satellite 
data,  pseudo-range  estiaates,  receiver  paraaeters,  and  position 
solution  estiaates  in  order  to  deteraLne  compliance  with  FAA 
navigation  requirements. 

c.  Microprocessor  Based  Design.  The  receiver  makes  extensive  use  of 
microprocessor  technology  to  synthesise  receiver  loops,  aanage 
receiver  operations  and  coapute  own  position. 

d.  Intelligent  Control  and  Display  Unit.  The  control  and  display 
uUt  is  managed  by  a  separate  Z&0  microprocessor  in  order  to 

allow  developmental  flexiUlity  in  the  pilot  display  Interface  design 
without  affecting  the  receiver  design. 

e.  Operational  Compatibility.  The  design  Incorporates  aodes  of 
operation  consistent  with  current  VOR  enroute  teralnal  and 
non-precision  approach  procedures.  In  addition,  it  includes  direct 
routing  area  navigation  (RNAV)  nodes. 

3.1.1  Systea  Tlalng 

The  typical  sequence  of  events  following  power-on  is  illustrated  in 
Fig.  3-2.  taring  a  brief  initialisation  period,  built-in  diagnostics  test 
microprocessor  and  receiver  functions,  and  the  age  of  the  stored  alaanac  are 
checked.  The  position  software  uses  the  last  stored  position  and  alaanac  data 
to  select  four  space  vehicles  (SV's)  for  initial  acquisition.  The  first  SV  is 
then  assigned  to  both  channels,  one  searching  the  1023  chip  code  in  a  1500  Hs 
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bin  centered  on  the  expected  frequency,  and  the  second  searching  in  an 
adjacent  1500  Hz  bin.  As  shown  in  the  exaaple  of  Fig.  3-2,  Channel  2  acquired 
SV-1  first  and  proceeded  to  deaodulate  the  1500  bit  navigation  message. 
Concurrently,  the  receiver  aanageaent  software  directed  Channel  1  to  search 
for  SV-2. 

When  four  SV's  are  acquired,  Channel  1  is  then  placed  in  the  transition 
node.  On  entry  to  the  transition  node,  the  satellites  have  ±  17  chip  code 
phase  and  ±  750  Hz  frequency  uncertainties.  After  several  four-second 
transition  cycles,  the  code  and  frequency  uncertainties  narrow  to  ±  0.16  chips 
and  ±  200  Hz.  Channel  1  is  then  assigned  to  the  navigation  node  where  SVs  1-4 
are  sequentially  accessed  at  220  msec  per  SV.  The  navigation  node  operates  on 
a  ten  slot  cycle  lasting  2.2  seconds.  As  additional  SVs  are  acquired  by 
Channel  2,  they  are  immediately  transferred  to  a  slot  in  the  Channel  1 
navigation  mode  cycle. 

When  six  SVs  have  been  acquired  or  six  minutes  has  passed,  a  smoothed 
position  estimate,  based  on  a  batch  processed  least-squares  linearization 
algorithm,  is  transferred  to  the  pilot  display.  It  is  Important  to  attenpt  to 
have  acquired  at  least  six  SVs  prior  to  activation  of  the  pilot  display  in 
order  to  have  two  backup  SVs  during  aircraft  takeoff  maneuvers.  Typically  six 
or  more  SVs  will  have  been  acquired  within  six  ninutes. 

3.1.2  Alternative  Startup  Modes 

The  startup  sequence  discussed  in  the  previous  section  assumed  that  the 
current  tine  (+  1  minute)  and  location  (+  3  miles)  are  available  (i.e.,  were 
recorded  prior  to  the  previous  power  shutdown)  and  that  a  stored  almanac 
exists  whose  age  is  less  than  30  days.  Should  own  position  be  unknown,  or  the 
battery  supporting  the  built-in  calendar  clock  be  low,  one  of  the  alternative 
startup  nodes  shown  in  Table  3-1  will  be  automatically  selected. 

Should  the  almanac  be  old  (age  >  30  days),  an  additional  12.5  minutes  is 
necessary  to  acquire  a  new  almanac  from  any  of  the  available  SVs. 

3.2  Antenna 

An  Important  aspect  of  the  low-cost  GPS  receiver  development  is  the 
definition  of  the  required  gain  versus  elevation  angle  characteristic  for  the 
aircraft  antenna.  Results  of  a  study  to  determine  this  characteristic  are 
summarized  below. 

3.2.1  General  Considerations 

It  is  convenient  to  consider  the  above-horizon  and  below-horizon  portion 
of  the  GPS  antenna  gain  characteristic  separately.  The  above  horizon 
characteristic  depends  on  the  satellite  slant  ranges  at  low  and  high  elevation 
angles.  Since  the  free  space  path  loss  is  2  dB  greater  at  5*  elevation  than 
at  zenith,  the  antenna  gain  should  be  2  dB  greater  at  low  elevations  than  at 
high  elevations.  The  boundary  between  these  two  regines  is  taken  as  30°, 
where  the  free  space  path  loss  is  1  dB  greater  than  at  5°. 


TABLE  3-1 


ALTERNATIVE  START-UP  MODES 


ASSUMPTIONS 


CHARACTERISTICS 


Cold  •  Full  range  and  L.O.  •  Prestored  almanac 

Start  frequency  uncertainty  (may  be  old) 

•  User  position  •  Stationary  user 

uncertainty 


•  Trial-and -error  SV 
visibility  determi¬ 
nation 

•  Hide  frequency 
searches 

•  Full  code  searches 

•  Constellation 
review  and  2nd 
acquisition  pass 
following  first  fix 

•  TTFF  -  10-36  mins 


“Time 

only" 

starts 


•  Know  time  to  i  1  ml n  •  Prestored  almanac 

•  Know  coarse  position  •  Stationary  user 
<±  6.75°) 

•  Full  L.O.  frequency 
uncertainty 


•  Predicted  SV  visi¬ 
bility 

•  Wide  1st  SV  frequency 
search  -  narrowed 
subsequently 

•  Full  code  searches 

•  TTFF  -  4-7.5  mins 


Time -and  •  Know  time  to  1  1  mLn  •  Prestored  almanac 
position 

wideband  •  Know  position  to  i  T  •  User  motion: 
start  200  lcnot  velocity 

•  Full  L.O.  freq  0.5g  acceleration 

uncertainty  30  deg  banks 

(3xl0-6) 


•  Predicted  SV  visibility 

•  Wide  1st  SV  frequency 
search  -  narrowed 
subsequently 

•  Full  code  searches 

•  TTFF  -  3-4  ndns 


Time-and  •  Know  time  to  ±  1  min 
position 

narrow-  •  Know  position  to  i  3 

band  miles 

start 


•  Know  L.O.  frequency 
to  ±  5  x  10“7 


•  Predicted  SV  visibility 

•  Narrow  frequency 
searches 

•  Full  code  1st  SV 
search  -  limited  sub¬ 
sequently 

s  TTFF  -  3-4  mins 


The  below-horizon  characteristic  depends  upon  the  effect  of  aultipath  on 
the  GPS  receiver.  Two  types  of  effect  result  from  aultipath.  The  first 
effect  is  that  deep  fading  of  the  GPS  signal  can  result  from  the  destructive 
Interference  of  the  direct  and  aultipath  signals.  The  oagnitude  of  this 
effect  is  largely  dependent  on  the  GPS  signal  fade  margin.  In  general,  a  deep 
fade  will  mean  either  a  temporary  loss  of  the  signal  or  a  delay  in  signal 
acquisition. 

The  second  effect  of  ailtipath  on  the  receiver  is  to  cause  an  error  in 
the  delay-lock-loop  pseudo-range  tracking.  The  nature  of  this  error  depends 
upon  whether  the  aultipath  delay  is  greater  or  less  than  1.5  C/A  code  chips. 
The  dependence  of  the  multipath  delay  upon  elevation  angle  and  aircraft 
altitude  is  shown  in  Fig.  3-3.  As  seen  in  the  figure,  the  multipath  delay 
always  exceeds  1.5  code  chips  for  altitudes  greater  than  8600  feet  and 
elevation  angles  greater  than  5°.  In  general,  tracking  errors  due  to 
aultipath  delays  greater  than  1.5  chips  can  be  reoioved  by  software  tests, 
whereas  errors  due  to  delays  less  than  1.5  chips  cause  tracking  bias  errors 
which  cannot  be  removed.  These  effects  are  discussed  more  fully  in  the 
following  paragraphs. 

3.2.2  Multipath  Delay  Greater  Than  1.5  Micro seconds 

If  the  multipath  delay  is  greater  than  1.5  nsec  (1.5  C/A  code  chips),  the 
delay-lock-loop  (DLL)  discriminator  characteristic  can  develop  a  second  stable 
operating  point,  as  shown  in  Fig.  3-4.  In  this  case,  the  DLL  nay  track  the 
pseudo-range  of  the  multipath  signal  rather  than  the  direct  signal,  resulting 
in  a  large  (>  1500  ft)  pseudo-range  error.  However,  since  the  receiver  is 
sequential,  reacquiring  each  satellite  once  per  second,  this  error  is  likely 
to  be  manifested  as  an  occasional  false  lock  with  accompanying  large  range 
change.  Thus,  the  receiver  control  software  can  eliminate  this  type  of  error 
by  tracking  the  pseudo-range  and  rejecting  any  unrealistically  large  range 
changes. 

3.2.3  Multipath  Delay  Less  Than  1.5  Microseconds 

In  the  case  where  the  ailtipath  delay  is  less  than  1.5  chips  the  effect 
on  the  DLL  performance  is  more  serious.  Figure  3-5  shows  that  the  effect  of 
aultipath  in  this  case  is  not  to  create  a  second  stable  operating  point,  but 
rather  to  Introduce  a  bias  into  the  original  operating  point.  This  bias, 
furthermore,  may  not  be  dlscernable  by  the  position  software.  The  magnitude 
of  this  error  depends  upon  both  the  magnitude  of  the  ailtipath  relative  to  the 
direct  signal  and  magnitude  of  the  aultipath  delay. 


Hageraan  (Ref.  9)  has  studied  the  effect  of  aultlpath  on  both  the 
coherent  and  noncoherent  Delay  Lock  Loop.  Two  regions  of  DLL  operation  In  the 
presence  of  aultlpath  have  been  defined.  In  Region  1,  the  DLL  Is  tracking  the 
direct  signal,  whereas  In  Region  II  the  DLL  Is  tracking  the  aultlpath  signal. 
The  results  of  Bageroan's  analysis  for  the  noncoherent  DLL  are  shown  In 
Figs.  3-6  and  3-7.  (Note:  Hageraan's  results  are  In  teras  of  the  P-code  and 
thus  oust  be  aultlplled  by  10  for  the  C/A  code).  Although  the  Region  II 
errors  are  auch  larger  than  the  Region  I  errors,  they  are  less  probable, 
especially  in  the  case  of  a  non-coherent  delay  lock  loop.  The  worst-case 
error  therefore  can  be  taken  froa  the  results  for  Region  I. 

It  Is  seen  froa  these  results  that  the  worst-case  error  depends  on  the 
relative  aultlpath  aaplitude.  If  the  aultlpath  aaplitude  can  be  Halted  to 
0.2,  then  the  expected  (aean)  tracking  error  is  10  feet  and  the  ras  error  Is 
50  feet.  The  aultlpath  aaplitude  depends  In  turn  on  the  antenna  gain  vs. 
elevation  characteristic  and  the  reflection  coefficient  of  the  surface. 
According  to  Figure  3-8  (Ref.  10)  the  worst  case  reflection  coefficient  at  5* 
is  0.7.  If  the  antenna  gain  is  0  dBIC  at  5*,  then  the  required  antenna  gain 
at  -5*  can  be  calculated  froa: 

Antenna  Gain  (4-5*) 

Relative  aultlpath  - - -— - x  Reflection  coefficient  at  5* 

aaplitude  at  5*  Antenna  Gain  (-5*) 

If  it  is  assuaed  that  the  antenna  gain  at  5*  is  at  least  0  dBIC,  then  the 
required  antenna  gain  at  -5“  should  be  no  greater  than  -5.5  dB. 

The  results  of  the  foregoing  discussion  are  suaaarixed  In  Table  3-2. 


TABLE  3-2 

Initial  Antenna  Gain  vs.  Elevation  Angle  Requlreaent 

Elevation  Angle  Antenna  Gain 

30*  to  90*  >  -2  dBIC 

5*  to  30*  >  0  dBIC 

-90*  to  -5*  <  -5.5  dBIC 


This  characteristic  Is  Illustrated  In  Fig.  3-9.  Also  shown  Is  the 
specification  developed  by  General  Dynaalcs  (Ref.  11).  The  single  difference 
between  the  two  characteristics  Is  that  General  Dynaalcs  specifies  <  -7  dBIC 
below  10*,  perhaps  due  to  the  difficulty  of  achieving  -5.5  dBIC  at  5*.  The 
receiver  perforaance  characteristics,  however,  will  be  essentially  equivalent 
for  both  specifications. 


2dbic 


FIG.  3-9.  INITIAL  ANTENNA  GAIN  VS  ELEVATION  ANGLE  REQUIREMENT 


However,  the  way  in  which  groups  of  pseudo-range  measurements  (up  to  10 
measurements,  from  as  many  satellites)  are  batch  processed  In  the  position 
computation  aakes  the  gain  characteristic  for  mltipath  suppression  less 
stringent.  The  batch  least-squares  position  estimation  method  reduces  the 
sensitivity  of  the  position  solution  to  an  Individual  aultlpath-corrupted 
pseudo-range  measurement.  With  least-squares  fit  processing  it  was  determined 
that  It  la  possible  to  redesign  for  a  random  error  of  120  feet,  as  long  as  the 
associated  bias  error  was  limited  to  50  feet,  worst  case.  The  below-horlson 
antenna  gain  requirements  were  recalculated  to  meet  this  revised 
multlpath-lnduced  range  error  bias  limit  assuming  a  ground  bounce  reflection 
coefficient  of  0.7,  maximum. 

The  revised  specifications  are  depicted  in  Table  3-3  and  in  Fig.  3-10. 

Note  that  the  below-horlson  gain  for  elevation  angles  in  the  range  minus  5  to 
minus  30  degrees  need  only  be  2  dB  less  than  the  gain  over  the  above-horlson  5 
to  30  degree  sector. 

TABLE  3-3.  REVISED  GPS  ANTENNA  SPECIFICATION 


Elevation  Antenna 

Angle,  (u)  Gain,  G(a) 


30* 

to 

90* 

> 

-2  dBIC 

5* 

to 

30° 

> 

0  dBIC 

-5* 

to 

5* 

< 

0  dBIC 

-30* 

to 

-5* 

< 

-2  dBIC 

-90* 

to 

-30* 

< 

-4  dBIC 

The  measured  performance  of  an  antenna  (Fig.  3-11)  purchased  from  Chu 
Associates  using  the  revised  specifications  is  reported  In  Section  4.2.1. 

3.3  Preamplifier 

The  receiver  preamplifier  specification  was  based  on  considerations  of 
noise  figure  and  out-of-band  rejection.  A  nominal  link  budget,  based  on 
reasonable  fade  margins,  resulted  in  the  need  for  a  4  dB  receiver  noise 
figure.  Since  3  dB  must  be  allocated  to  the  preamplifier  and  1  dB  to  losses 
associated  with  filters  and  connectors.  It  is  necessary  that  the  preamplifier 
be  located  at  the  antenna.  The  dual-channel  receiver  developed  by  STI  further 
assumed  a  nominal  gain  of  50  dB  with  10  dB  cable  loas  between  the  amplifier 
and  receiver. 

A  commercial  amplifier,  an  Avantek  AM-1664M,  was  purchased  with  the 
following  specifications: 

fQ  1575.42  MHx 

Bandwidth  20  MHs  (3  dB) 

52  dB 


Gain 

Noise  Figure 


2  dB 


oiap  i 
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3.4  Dual  Channel  Receiver  Hardware 


3.4.1  Architecture  and  Design  Features 

The  receiver,  developed  by  Stanford  Telecommunications,  Inc.  for  Lincoln 
Laboratory,  consists  of  two  Identical  channels  fed  by  a  common  down-converter 
and  timer/ synthesizer  as  shown  in  Fig.  3-12.  Each  channel  performs  digital 
functions  in  Z8000  microprocessors.  The  functional  schematic  of  Fig.  3-13 
provides  further  detail  for  one  of  the  two  channels  with  analog  (hardware)  and 
digital  (software)  functions  to  the  left  and  right  of  the  dashed  vertical 
line.  Software  switches  SHI  and  SW2  are  programmed  so  that  during  acquisition 
SHI  is  in  the  “selected  frequency"  position  and  SH2  in  “code  search".  During 
the  navigation  mode  a  steady  state  condition  will  occur  in  which  SHI  is 
programmed  to  “AFC"  and  SH2  to  the  "DLL"  Position.  The  allocation  of  receiver 
functions  between  the  Z8000  microcomputers  and  the  position  software  in  the 
LSI  11/23  is  shown  in  Table  3-4. 

The  IF  noise  bandwidth,  BIF,  is  determined  by  the  integrate  and  dump  (I 
and  D)  circuits  and  by  the  digital  accumulators  according  to  the  expression 


®IF  - 

2  T  M 

where  T  is  the  sampling  period  and  M  the  number  of  sables  accumulated.  The 
output  of  the  low-pass  filter  is  sampled  at  a  2  KHz  rate.  Additional 
filtering  is  provided  by  the  post-detection  Integrator  (PDI). 

The  code  select  circuitry  can  operate  in  the  punctual  mode  (P)  for  code 
search  and  data  extraction,  or  in  the  early-late  dither  mode  to  control  the 
delay  locked  loop  (DLL).  Characteristics  of  the  search  detector,  and  the  AFC, 
code  and  AQC  loops  are  shown  in  Table  3-5. 

Figure  3-14  shows  the  AFC  lock  detector  and  the  carrier  and  noise  ratio 
estimator  embodied  in  each  STI  receiver  channel.  Characteristics  of  these 
circuits  are  given  in  Table  3-6. 

3.4.2  Receiver  States 

When  a  receiver  channel  is  provided  with  commands  and  initialization  data 
by  receiver  management  software  in  the  LSI  11/23  it  responds  by  executing  one 
of  the  following  Z8000  programs: 

•  initialization 

•  acquisition 

•  transition 

•  navigation 

•  data  demodulation 
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FIG.  3-10.  FUNCTIONAL  SCHEMATIC  OF  RECEIVER 


TABLE  3-4 

RECEIVER  CHANNEL  -  POSITION  SOFTWARE  FUNCTIONAL  ALLOCATION 


REQUIREMENT 

Z8000 

RECEIVER  CHANNEL  SOFTWARE 

LSI  11/23 

POSITION  SOFTWARE 

Acquistion 

Sequential  detection 

Pull-In  and  track 

Data  synchronization 

SV  selection 

Acquisition  prepositioning 
Search  control 

Acquisition  sequence  control 

Sequential 

Track 

Limited  sequential  detection 
Pull-in  and  track 

Dwell  sequencing 

Inter-dwell  prepositioning 

SV  constellation  revision 

New  SV  prepositioning 

Sequence  monitoring  and 
control 

Paeudo-range 

Code  state  sampling 

Pseudo-range  construction 
Ambiguity  resolution 
Ionlspherlc  and  tropospheric 
compensation 

Data 

Demodulation 

Bit  demodulation 

Word  recovery 

Parity  checking 

Data  synchronization  (single 
channel) 

Data  unpacking 

Demodulation  monitoring  and 
control 

TABLE  3-5 


RECEIVER  LOOP  CHARACTERISTICS 


SEARCH  DETECTOR 

AFC  LOOP 

CODE  LOOP 

AGO 

Square  law  power 

Delayed  cross 

1  0.5  chips  early/ 

Maintains 

detector 

product  detector 
«  sin  2nAf  x 

late  power  detector 

constant  power 

Open  loop 

1st  order  loop 

Carrier  (AFC  or  PLL) 

Deals  with  up 

operation 

loop  aided 

to  10  satellite 
signals . 

AGO  sets  noise 

Variable  loop 

Variable  loop 

Variable  loop 

floor 

bandwidth 

bandwidth 

bandwidth 

Two  search  rates: 

Primary  loop 

Primary  loop 

Primary  loop 

50  chips /sec 

bandwidth  5  Hz 

bandwidth  5  Hz 

bandwidth  1  Hz 

83.3  chips/sec 

and  10  Hz 

and  10  Hz 

(for  1st  satellite 
acquisition  only). 
10  Hz  otherwise 

Multiple  freq. 

Operates  with 

cell  search 

early/late  code 

Supplies  power 
reference  to  lock 
detectors 

Upon  detection  does 
false  alarm  check 

TABLE  3-6 


RECEIVER  AFC  LOCK  DETECTOR  AND  C/N0  ESTIMATOR  CHARACTERISTICS 


AFC  Lock 

Detector 

Detector  output 
■  cos  2*Afx 

Measures  AGC 
normalized  power 

Uses  narrow  band 

IF  (Bif-50  Hz)  for 
narrow  frequency 
discrimination 

Output  monotonic 
with  C/N0 

AGC  compensated  to 
be  insensitive  to  noise 
and  gain  variations 

Estimate  for  each 
satellite 

Also  serves  as  code 
sync  detector 

Filtered  over 
multiple  intervals 

Detection  bandwidth 
adjustable  from  5  Hz  up 

Detection  bandwidth 
5  Hz 

33 


ANALOG  .  DIGITAL 


The  Initialization  routine  establishes  operating  paraaeters  and  checks 
statue.  The  acquisition  routine  is  illustrated  in  Fig.  3-15.  Note  that  the 
IF  bandwidth,  Bjp,  is  changed  several  times  during  the  acquisition  sequence, 
which  illustrates  the  significant  advantage  of  software  realisations  of 
receiver  detection  algorithas. 

Figures  3-16  and  3-17  show  the  transition  and  navigation  state  diagraas. 

Figure  3-18  shows  the  detailed  events  that  occur  during  each  220 
Billisecond  dwell  Internal  in  navigation  aode.  The  prepositioning  data  will 
have  a  ±  0.16  chip  uncertainty  for  an  SV  entering  the  dwell  froa  phase 
transition.  Should  lock  fail  to  occur  twice,  the  receiver  will  autoaatically 
coaance  all  chip  search.  If  two  additional  failures  occur,  a  1  2  chip 
search  is  Initiated.  After  a  total  of  6  failures,  the  receiver  aanageaent 
software  (in  the  LSI  11/23)  will  coaaand  a  broader  search  by  relocation  of  the 
search  window.  After  several  subsequent  unsuccessful  atteapts  the  SV  will  be 
declared  to  be  potentially  faded  and  less  aggressive  techniques  will  be 
applied  to  reacquire  it. 

During  normal  operation  new  SVs  will  become  visible.  Figure  3-19  shows 
the  state  diagraas  for  the  acquisition  and  data  demodulation  of  a  new  SV. 

The  receiver  hardware  was  physically  asseabled  in  an  ATR  Chassis  as  shown 
in  Figs.  3-20  through  3-22.  The  receiver  modules  are  aounted  in  a 
multibus -type  chassis,  and  consist  of  a  mixture  of  off-the-shelf  commercial 
and  custom-built  boards.  It  should  be  noted  that  the  actual  size  of  a 
comae rcial  equivalent  will  be  significantly  (-8-10  tiaes)  sasller  because: 

a)  An  extra  board  to  provide  the  9K  byte  ROM  aeaory  for  each  Z8000  was 
Included  (the  Z8000  boards  provided  only  8K  of  local  ROM) . 

b)  LSI  devices  are  currently  being  developed  for  the  C/A  Coder  and  NCO 
functions  and, 

c)  much  of  the  board  area  is  not  populated. 


TIME 

BUDGET  REMARKS 


1  SEC 

(FIRST  SATELLITE 
ONLY) 


MODE,  SATELLITE 
IMPENDENT 

20.66  sec  xlCELLS 
+  0.60  -  1.12  sec. 

18  ns. 


^  -  1  Hs  (AGC  LOOP  BANDWIDTH) 

AGC  LOOP  SETTLES 
B^  -  1000  Hs 

N  CELLS  AS  REQUIRED 
PLUS  FALSE  ALARM  CHECK 

Bir  -  500  Hs  (IF  BANDWIDTH) 

3  CELLS,  1  CODE  POSITION 


60  ms. 


140  ms. 


Bxr  -  1000  Hs 

*LC  "  4  Hz,  -  10  Hz 

Bjy  -  500  Hz 

Bj*  “  2  Hz,  Bla  “  10  Hz 


220  ms 


2  sec 


7.22  -  13.22  sec 


Bj,  -  500  Hz 

®LC  “  2  ®LA  "  10  112 

Bjy  -  500  Hz 

®LC  “  2  Hz,  Bla  -  10  Hz 

■jf  "  50 

®LC  “  2  Hz,  Bla  -  10  Hz 


10.8  -  22.8  sec.  B^  -  50  Hs 

BLC  “  2  Hz,  -  10  Hz 

SEARCH  RATE:  50  CHIPS  PER  SECOND 


m.  *-16.  Receiver  Acquisition  State  Diagri 


enter  transition 


MODE,  SATELLITE 
DEPENDENT 

680-740  ms . 


20  ms. 


60  ms. 


140  ms. 


BIf  -1000  Hz 

+  15  CHIP  SEARCH 

BIp  -  500  Hz 
2  TRIES 


BIf  -  1000  He 

BLC  ‘  4  Hz*  BLA 

Bjp  -  200  Hz 

BLC  "  2  bla 


SEARCH  RATE:  50  CHIPS  PER  SECOND 
*LC  "  CODE  LOOP  BANDWIDTH 
Bn  -  AFC  LOOP  BANDWIDTH 


FKL  *-i*.  Rscslvsr  Transition  Stata  Diagram 


SLEW  CODERS, 
WAIT  FOR  BIT  EDGE 


L_ 

AFC,  CODE  LOOP 
PULL-IN 

~r 

AFC  LOCK 
DETECTION 


(  EXIT  ) 


TIME 

BUDGET 

<20  aS 


80-100  aS 


Bjf  -  200  Hz 
BLC  *  BLA  “  10  H 


100  aS 


Bip  ■  200  Hz 
BLC  *  "  5 


Bl(,  •  CODE  LOOP  BANDWIDTH 
B^-  AFC  LOOP  BANDWIDTH 


no.  9-17.  Racaivar  Navigation  State  Dlagr. 


TIME 

BUDGET 

6.3  SEC 
(±150  CHIPS) 


1000  Hz 
PLUS 


FALSE  ALARM  CHECK 


18  ms. 


60  ms. 


140  ms. 


220  ms. 


2  sec. 


B__  -  500  Hr 

3  CELLS f  1  CODE  POSITION 
PLUS  FALSE  ALARM  CHECK 

Bif  «  1000  Hz 

BLC  “  4  Hz’  BLA  "  10  112 

BJp  -  500  Hz 

®LC  "  2  Hz*  BLA  "  10  Hz 

BIF  ■  500  82 

blc  -  2  Hz>  bla  -  10  Hz 

Bif  *  500  Hz 

®LC  "  2  HZ’  ®LA  "  10  Hz 


7.22  -  13.22  sec 


®IF  "  50  Hz 

®LC  ‘  2  Hz’  BLA  “  10  Hz 


10.8  -  22.8  sec. 


Blp  -  50  Hz 

blc  ■  2  »•  “u  ■  10 
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3. 5  Position  Software 


3.5.1  Overview 

The  GPS  T  and  E  Equipment  software  is  organized  into  four  functional 
areas  as  shown  in  Fig.  3-23.  These  functional  areas  are  described  in 
Table  3-7.  Section  3.4  described  receiver  functions  resident  in  the  Z8000 
nl crocoaputers -  The  Navigation,  Instrumentation  and  Control  and  Display  Unit 
software  designs  are  described  in  Sections  3. 6-3.8. 

This  section  describes  the  position  software  located  in  the  LSI-11/23 
Position  Processor.  The  position  software,  developed  by  Intermmtrlcs ,  Inc., 
is  organized  into  several  functional  areas  as  shown  in  Fig.  3-24.  The 
following  paragraphs  describe  each  module: 

3.5. 1.1  System  Mode  Control 

This  module  maintains  control  over  all  other  modules.  It  controls  the 
receiver  channel  modes,  prepositions  the  receiver  channels  in  frequency  and 
code  phase,  and  updates  the  SV  list  as  visibilities  change  with  time  and 
aircraft  altitude. 

3.5. 1.2  System  Time  Management 

The  system  time  management  module  establishes  and  maintains  system  time 
based  on  the  50  Hz  clock  sent  from  the  dual  channel  receiver.  After 
correction  to  GPS  time  following  acquisition,  it  maintains  time  for  data 
tagging  and  display  of  GMT  to  the  operator  or  pilot. 

3.5. 1.3  Satellite  Selection 


This  module  makes  use  of  a  stored  almanac,  stored  current  position,  and  a 
battery-operated  calendar  clock  to  select  four  SVs  for  initial  acquisition. 

It  also  will  implement  a  search  strategy  should  less  orbital  information  be 
available  or  the  clock  have  failed.  While  in  operation  it  also  maintains  a 
list  of  up  to  10  visible  SVs. 

3. 5. 1.4  Receiver  Management 

The  functional  task  breakdown  between  the  position  software  and  the  dual 
channel  receiver  was  shown  in  Table  3-4.  A  high  level  flow  chart  which 
describes  the  major  receiver  states  from  the  point  of  view  of  the  receiver 
firmware  is  provided  in  Fig.  3-25.  The  detailed  activities  of  this  module 
during  acquisition  were  shown  in  Fig.  3-2. 

3. 5. 1.5  GPS  Navigation  Data  Management 

This  module  receives  parity-checked  navigation  message  words  from  the 
dual  channel  receiver.  It  monitors  the  age  of  the  stored  SV  ephemerls  data, 
requests  updates,  monitors  data  collection  and  maintains  a  full  data  almanac 
for  all  SVs. 


FKL  3-23.  3Y8TEM  SOFTWARE 


TABLE  3-7 


SOFTWARE  FUNCTIONAL  AREAS 


•  Receiver  Firmware  (Z8000) 


Directs  detailed  acquisition  of  SVs,  maintains  receiver 
feedback  loops  (AFC,  DLL,  etc.),  computes  pseudo-range, 
demodulates  navigation  data  and  maintains  dwell  tracking.  Is 
under  the  control  of  Receiver  Software. 

•  Position  Software  (LSI-11/23) 

Selects  satellites,  controls  Z8000  receiver  activities, 
computes  position  estimates  in  latitude-longitude  coordinates  and 
monitors  performance. 

•  Navigation  and  Instrumentation  Software  (LSI-11/23) 

Accepts  position  estimates  from  the  receiver  management 
software  and  waypoint  definitions  from  the  cockpit  display. 
Computes  bearing,  distance  and  time  estimates,  controls  course 
deviation  Indicator  and  provides  data  to  the  control  and  display 
unit  software.  Collects,  distributes  and  records  data. 

•  Control  and  Display  Unit  Software  (Z80) 

Establishes  specific  formats  and  protocol  for  data  entry  and 
display  on  the  cockpit  display  and  control  unit. 


LSI  11/23  POSITION  PROCESSOR 


FIQ.  3-24.  POSITION  SOFTWARE  SY8TEM  BLOCK  DIAGRAM 


3. 5. 1.6  Measurement  Processing 


This  nodule  collects  neasurenents  fron  the  receiver  channels,  resolves 
neasurenent  ambiguities,  compensates  for  deterministic  error  effects  such  as 
tropospheric  delay,  and  computes  a  position  and  time  fix  every  2.2  seconds. 

3.5.1.7  Navigation  Tracker 

This  module  maintains  a  current  estimate  of  the  aircraft  position, 
smooths  position  fixes  using  an  alpha-beta  tracker  and  interpolates  the  data 
to  provide  1-second  updates  to  the  navigation  software. 

3.5.1.8  DAD  Interface 

This  module  provides  linkage  with  the  the  Data  Acquisition  Unit  (DAU). 

It  accepts  inputs  during  system  start  up  (position  and  time  if  required),  and 
transfers  smoothed  navigation  estimates  (in  latitude-longitude  coordinates) 
and  aircraft  velocity  estimates  to  the  navigation  software  in  the  DAD.  The 
DAD  Interface  also  accepts  receiver  channel  configuration  commands  from  the 
DAO,  and  provides  receiver  data  to  the  DAD  for  recording  on  magnetic  tape. 

The  DAD  also  provides  a  limited  data  display  capability  while  in  flight.  The 
interface  to  the  DAD  is  a  DMA  hardware  Interface  since  the  DAU  resides  in  a 
separate  LSI-1 1/23. 

The  implementation  of  the  receiver  software  is  under  the  Digital 
Equipment  Corporation  RSX-11S  operating  system.  The  software  was  written  in 
RATFOR,  a  structured  FORTRAN  language.  A  top-down  design  approach  was  used  to 
achieve  developmental  flexibility  and  modular,  well  documented,  maintainable 
software. 

3.5.2  Acquisition  Strategy 

The  basic  acquisition  strategy  was  described  in  Section  3.1.  This 
section  provides  additional  detail  regarding  satellite  selection,  receiver 
channel  control  and  preparation  for  entry  into  the  transition  mode. 

While  provisions  for  the  cold  start  mode  were  Incorporated  in  the  design, 
the  typical  situation  as  described  here  will  be  a  narrow  or  wide  band  start  as 
defined  in  Table  3-1.  Both  assume  that  a  battery  operated  clock  provides  GMT 
to  *  1  minute,  and  that  the  present  location  and  an  almanac  less  than  30  days 
old  are  available  from  a  non-volatile  memory. 

A  satellite  selection  algorithm  is  used  in  the  acquisition  mode  in  order 
to  a)  select  a  reasonable  set  of  up  to  6  satellites  to  provide  good  geometry 
and  shadowing  protection  during  departure  if  one  or  two  satellites  fade  during 
turns,  and  b)  present  an  acceptable  fix  to  the  CDD  within  6  minutes.  Later, 
in  the  navigation  mode,  selection  is  unnecessary  as  all  satellites  in  view  are 
tracked. 


Initially,  a  set  of  four  satellites  are  selected  based  on  naxinising  the 
sua  of  the  azinuth  and  elevation  differences  between  all  pairs,  and  providing 
at  least  30*  elevation  angle  to  all  satellites,  if  possible.  Additional 
satellites,  if  available,  are  selected  for  subsequent  acquisition  in  order  to 
enter  the  navigation  node  with  up  to  6  in  track.  Eventually  all  satellites  in 
view  are  acquired  and  placed  in  the  navigation  node. 

The  tine  to  acquire  the  first  satellite  was  shown  by  Interaetrics  to  have 
a  snaller  variance  if  both  channels  searched  in  parallel  for  the  first 
satellite  than  if  each  channel  searched  independently  for  different 
satellites*  Therefore,  initial  acquisition  is  done  by  both  channels  searching 
cooperatively  for  a  single  satellite.  When  the  first  satellite  is  acquired, 
the  local  clock  is  corrected  to  a  value  within  20  Billiseconds  of  GPS  tine  and 
the  reference  oscillator  frequency  error  is  logged  to  reduce  the  frequency 
search  during  the  remaining  satellite  acquisition. 

The  acquisition  strategy  includes  reacquiring  all  previously  acquired 
satellites  following  each  new  acquisition  (a  detail  not  shown  in  Fig.  3-2) 
because  of  the  need  for  recent  code  phase  and  doppler  estinates  in  order  to 
preposition  the  receiver  channel  for  the  transition  node.  If  acquisition  is 
conpleted  but  only  three  satellites  are  located,  the  barometric  altitude  is 
substituted  and  the  initial  position  fix  is  computed. 

3.5.3  Transition  Strategy 

When  four  satellites  have  been  successfully  acquired  (or  if  necessary, 
three  plus  barometric  altitude),  the  position  software  connands  one  of  the 
receiver  channels  into  the  transition  node.  In  this  node,  each  satellite  ia 
acquired  and  tracked  for  1  second.  At  the  end  of  each  four  second  schedule  a 
fix  is  attempted.  When  two  successive  fixes  are  successful,  the  receiver 
software  coaputes  prepositioning  data  in  preparation  for  the  navigation  node. 

3.5.4  Navigation  Strategy 

When  the  position  software  determines  that  two  successive  fixes  have  been 
successful,  the  transition  node  is  terminated  and  the  receiver  channel  placed 
in  the  navigation  node.  Once  started  the  receiver  channel  Z8000  is  able  to 
sustain  the  0.22  second  per  satellite  sequential  detection  of  up  to  ten 
satellites  without  further  assistance  fron  the  position  software.  The 
position  software  does,  however,  have  to  synchronously  request  and  accept 
pseudo-range,  code  phase,  and  status  data  fron  the  navigation  node  channel 
every  0.22  seconds. 

Once  one  receiver  channel  is  placed  in  the  navigation  node,  the  second 
channel  is  used  to  acquire  and  track  additional  satellites  in  order  to  place 
then  in  a  navigation  node  slot,  and  to  update  ephenerls  and  almanac  data.  As 
new  satellites  rise,  a  satellite  visibility  algorithm  (run  every  3  minutes) 
enables  SV  acquisitions  by  the  second  channel  and  allows  new  SVs  to  be  placed 
in  navigation  slots  in  the  first  channel. 


The  ten  navigation  slots  are  all  filled  if  less  than  ten  satellites  are 
in  view,  by  repeating  satellites  in  the  extra  slots  as  necessary.  This 
simplifies  several  algorithms  and  provides  additional  margin  for  satellites 
that  are  fading  but  are  assigned  to  more  than  one  slot. 

A  pseudo-range  measurement  is  declared  valid  by  the  receiver  channel  if 
AFC  lock  occurred  during  the  last  dwell  interval,  and  it  is  within  1000  feet 
of  a  predicted  value  based  on  a  predicted  position.  When  a  valid  measurement 
is  not  received  for  a  particular  satellite,  the  algorithm  in  Fig.  3-26  is 
invoked. 

Valid  pseudo-range  measurements  for  each  2.2  second  scheduling  cycle  are 
propagated  to  a  common  solution  time  using  measured  doppler  and  presented  to 
the  position  estimation  algorithm  described  in  Section  3.5.5. 

3.5.5  Position  Estimation 

A  batch-processed  linearization  (BPL)  position  fixing  algorithm  was 
selected  following  a  study  of  alternative  techniques.  The  principal 
advantages  of  the  BPL  algorithm  are  that  a)  it  uses  all  available  pseudo-range 
measurements  in  every  position  fix,  b)  it  is  a  minimum  least-squares  error 
estimate,  and  c)  it  can  be  easily  Implemented.  A  description  of  the  position 
fixing  algorithm  is  provided  in  Appendix  A. 

A  separate  analysis  of  tracking  filters  was  conducted  which  compared 
alpha-beta,  least-squares  second  order,  and  Kalman  filters.  The  functions  of 
the  filter  are  to  a)  accept  new  position  fixes  every  2.2  seconds,  b)  compute 
new  smoothed  position  estimates  at  a  1  Hz  rate,  and  c)  coast  the  output  for  up 
to  8  seconds*  when  new  position  estimates  do  not  occur.  Should  no  estimates 
be  available  for  8  seconds,  the  NAV  flag  is  shown  to  the  pilot  indicating 
unacceptable  performance. 

The  emphasis  on  low  cost  and  limited  aircraft  dynamics  (200  KTS,  0.5g) 
led  to  the  selection  of  a  fixed  gain  alpha-beta  filter  which  had,  in 
simulation,  acceptable  dynamic  performance  and  the  lowest  computational 
requirements.  The  relationship  between  the  position  fixing  algorithm  and  the 
navigation  tracker  is  Illustrated  in  Fig.  3-27.  A  more  detailed  discussion  of 
the  navigation  tracker  is  provided  in  Appendix  B. 

3.5.6  Performance  Monitoring 

Performance  monitoring  features  are  provided  in  each  receiver  channel  and 
in  the  position  software.  Table  3-8  lists  the  various  parameters  that  are 
monitored  in  each  receiver  channel. 


*  Selected  to  limit  the  horizontal  position  error  to  500  feet 


AGO  gain  out  of  range 
C/A  coder 

Loss  of  phase  lock  on  synthesizer  interrupt 

20-msec  system  clock  not  present 

NCO  frequency  word  greater  than  maximum  limit 

Correlator  power  not  in  range 

Correlator  2-KHz  Interrupt  not  present 

CPU  fault 


The  basic  performance  monitoring  algorithm  in  the  position  software  is 
shown  in  Fig.  3-28  and  3-29. 

The  position  software  also  tests  the  navigation  data  for  inconsistent 
values  such  as  sudden  change  in  clock  bias.  Slow  drifts  are.  however,  more 
likely  and  may  be  difficult  to  detect  if  operating  with  a  mini mun  number  of 
satellites  and  high  GDOP. 

3.5.7  Fail-Soft  Techniques 

Several  fail-soft  techniques  were  developed.  If  during  any  mode  one  of 
the  two  receiver  channels  fails,  a  design  was  developed  in  which  the  remaining 
healthy  channel  provides  reduced  rate  navigation  service  to  the  pilot 
displays.  The  single  channel  case  operates  in  one  of  two  modes: 

a)  A  mode  identical  to  the  high  rate  navigation  mode  of  the  dual  channel 
case. 

b)  A  nixed  mode,  illustrated  in  figure  3-30,  which  allows  data 
demodulation  to  occur  for  one  satellite  and  six  satellites  to  be 
tracked  on  a  2.4  second  cycle. 

The  latter  mode  requires  a  satellite  selection  algorithm  to  select  the 
six  satellites  for  the  navigation  dwells.  The  mixed  mode  has  been  fully 
implemented  and  validated  in  the  receiver  channels,  and  designed  but  not 
lnplemented  in  the  position  software. 

As  indicated  in  the  description  of  the  position  fixing  algorithm.  Section 
3.5.5,  the  position  software  monitors  HOOP  and  the  number  of  satellites  in 
track.  As  a  result  it  takes  the  actions  summarized  in  Table  3-9. 


Another  fail-soft  technique  addressed  the  problem  of  rapid  recovery 
following  a  momentary  power  outage.  An  algorithm  was  designed  which  restored 
navigation  service  within  one  minute  of  the  end  of  a  power  outage  lasting  less 
than  30  seconds.  The  technique  requires  that  receiver  data  be  remembered  and 
that  a  battery  backup  be  provided  for  the  reference  oscillator  (to  sustain  the 
oscillator  and  system  clock  for  30  seconds).  This  technique  was  not 
lflg>lemented  due  to  other  higher  priority  work  and  because  it  appeared  to  be 
straightforward . 

3.6  navigation  Software 

The  navigation  software  is  organized  as  shown  in  Fig.  3-31.  It  receives 
smoothed  position  estimates  from  the  position  software  at  a  once-per- second 
rate.  The  Area  Navigation  Computation  (RNAV)  module  uses  the  position 
estimates  to  compute  navigation  data  for  display  to  the  pilot.  As  shown  in 
Fig.  3-32,  this  includes  dlstance-to-waypolnt  (DIST),  bearing  to  waypoint 
(BRG)  and  cross-track  deviation  (XTK) .  The  cross-track  deviation  is  displayed 
to  the  pilot  via  the  course  deviation  indicator  shown  in  Fig.  3-33;  the  other 
parameters  are  displayed  via  the  CDU  front  panel. 


NOTE:  REACQUISfTION  OF  3V  REJECTED 
FOR  BAD  PSEUDO  RANGE  18 
ATTEMPTED  AFTER  SOO  SECOND 
DELAY. 


FIG.  3-28.  POSITION  SOFTWARE  PERFORMANCE  MONITOR 
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FIG.  3-29.  POSITION  SOFTWARE  PERFORMANCE  MONITOR 
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FIG.  3-30.  Timing  Diagram  for  a  2.4  Second  Cycle  Single  Channel  System. 
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Declare  Failure; 
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3 
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altitude  (Pseudo-baro). 


3 
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Baromentric  Altitude 
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FIG.  3-33.  COUR8E  DEVIATION  INDICATOR/OMNI-BEARING  SELECTOR 


The  CDU  interface  (CDUI)  nodule  receives  navigation  data  fron  the  RNAV 
nodule  and  passes  the  data  to  the  CDU  for  front  panel  display.  CDUI  also 
receives  pilot  data  fron  the  CDU  consisting  of  such  data  as  startup  node,  CDI 
sensitivity,  and  active  waypoint  nunber.  The  RNAV  nodule  passes  the 
cross-track  deviation  to  the  CDI  and  reads  the  desired  course  fron  the  OBS. 

3.7  Control  and  Display  Unit 

The  function  of  the  Control  and  Display  Unit  (CDU)  is  to  interface  the 
pilot  to  the  test  and  evaluation  equipnent.  The  CDU  front  panel  layout  is 
shown  in  Fig.  3-34.  The  front  panel  consists  of  an  alphanumeric  display  with 
two  lines  of  16  characters  per  line,  a  32  position  keyboard  and  three  control 
switches.  The  alphanunerlc  displays  are  of  the  LED  lighted -segment  type  and 
are  0.25  Inches  high.  These  displays  are  bright,  easily  read  at  a  distance  of 
five  feet,  and  feature  a  wide  viewing  angle  (+  55°).  The  keypads  are  compact, 
off-the-shelf,  units  with  0.5  inch  button  spacing;  control  switches  are 
standard  rotary  types. 

CDU  functions  can  be  divided  into  two  groups:  flight  plan  data  entry  and 
navigation  data  display.  Flight  plan  entry  occurs  before  takeoff  while  the 
GPS  receiver  is  still  in  the  acquisition  mode  and  has  not  yet  obtained  an 
initial  position  fix.  The  pilot  uses  this  time  to  enter  the  desired  flight 
plan  via  the  CDU  keyboard,  with  appropriate  prompting  messages  from  the 
alphanumeric  display.  The  flight  plan  data  is  passed  to  the  Navigation 
Software  via  an  ARINC  429  Interface,  as  shown  in  Fig.  3-35.  A  sample  flight 
plan  for  a  trip  from  Atlantic  City,  NJ  to  Hanscom  Field,  HA  is  illustrated  in 
Fig.  3-36.  The  flight  plan  data  is  entered  in  the  form  of  waypoints  along  the 
planned  flight  path.  The  data  file  generated  by  this  data  entry  process  is 
termed  the  Stored  Waypoint  data  file.1  The  Stored  Waypoint  data  file  for  the 
sample  flight  plan  is  shown  in  Table  3-10. 

The  waypoints  along  the  flight  plan  can  be  of  several  types:  VDR, 
intersection,  latitude/longitude  or  artificial.  VDR  and  intersection 
waypoints  are  entered  by  three-and  five-letter  labels,  respectively. 

Artificial  waypoints  are  specified  in  terns  of  range  and  bearing  fron  a  VOR, 
intersection  or  latitude/longitude.  As  the  waypoints  are  entered,  they  are 
accumulated  in  the  Stored  Waypoint  data  file.  The  navigation  software 
determines  the  latitude  and  longitude  of  each  waypoint  using  look-up  tables 
and  calculations  as  required.  The  Stored  Waypoint  data  file  also  Includes  the 
courses  to  and  from  each  waypoint.  These  courses  can  be  entered  by  the  pilot 
or  computed  by  the  navigation  software. 

The  other  function  of  the  CDU  is  to  display  navigation  data  while  the 
aircraft  la  in  flight.  The  data  required  to  support  the  navigation  display  is 
termed  the  Active  Waypoint  data  file.  This  data  file  Is  maintained  by  the 
navigation  software  and  passed  to  the  CDU  via  the  ARINC  Interface.  The  Active 
Waypoint  data  file  must  be  updated  at  a  rate  such  that  the  CDU  navigation 
display  is  always  current  (l.e.,  every  second). 


1.  Due  to  schedule  constraints,  the  flight  plan  data  entry  function  was  not 
implemented  in  the  navigation  software.  Instead,  the  stored  waypoint  file  was 
predefined  in  the  navigation  sof tv  —  *-o  include  a  set  of  waypoints  for  the 
test  flight  area.  The  effect  was  the  sane  as  if  the  pilot  had  entered  the 
waypoints  manually. 
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CONTROL  AND  DISPLAY  UNIT  (CDU) 


DAU-CDU  INTERFACE 


WFT 

TYPE1 

IDBNT 

TO 

FROM 

LATITUDE 

LONGITUDE 

OFFSET 

0 

L/L 

ACY2 

- 
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39*27. 4'N 
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1 
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59 

Lookup 
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3 

VOR 
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56 
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CHOMR 

72 
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VOR 

ORW 

11 

30 

— 

M 

M 

6 

VOR 

POT 

30 

25 

•9 
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•S 

7 

AWP 

V4312 

25 

111 

Calculate 

Calculate 

GDM  111* 
16.0  nni 

8 

INT 

LOBBY 

111 

113 

Lookup 

Lookup 

None 

9 

NDB 

BE 

112 

112 

M 

*• 

M 

10 

L/L 

BED2 

112 

114 

42*28. 2'N 

71*17.4'W 

None 

11 

NDB 

SKR 

114 

294 

Lookup 

Lookup 

None 

Not*  1:  L/L  “  Latltude/Longitud* 

I NT  “  Intersection 

AWP  -  Artificial  Waypoint 

NDB  ■  Non-dlrectlonal  Beacon 


Note  2:  IDSNT  Optional  for  L/L  and  AWP  types 


The  contents  of  the  Active  Waypoint  data  file  are  shown  In  Table  3-11. 
This  data  Includes  such  Information  as  the  range  and  bearing  to  the  active 
waypoint,  i.e.,  the  waypoint  the  aircraft  is  currently  heading  to  or  from. 

The  active  waypoints  can  be  selected  by  the  pilot  via  the  CDU  keyboard  or  can 
be  automatically  incremented  as  the  waypoints  along  the  flight  plan  are 
overflown.  The  specific  data  displayed  is  selected  by  keyboard  entries  as 
summarised  in  Table  3-12.  The  navigation  data  appears  in  the  "DISPLAY"  field 
of  the  CDU  display. 

A  typical  CDU  navigation  display  is  shown  in  Fig.  3-37.  The  display 
indicates  that  the  CDU  mode  is  "OPERATE",  the  active  waypoint  is  number  6 
(Putnam  VOR),  and  that  the  aircraft  is  19.2  nmi  away  from  the  waypoint  and 
traveling  toward  it.  Alternate  data  can  be  obtained  by  pushing  one  of  the 
keys  in  the  leftmost  three  columns.  For  example,  pushing  "TTG"  would  display 
the  time-to-go  to  the  active  waypoint.  The  active  waypoint  could  be  changed 
to  number  5  by  pushing  "WPT",  "5"  and  "ENT". 

The  uppermost  rotary  switch  selects  the  CDI  sensitivity  desired.  The 
three  settings  are  for  Enroute,  Terminal  and  Approach,  in  order  of  least  to 
most  sensitive.  The  sequencing  switch  controls  whether  or  not  the  active 
waypoint  will  increment  automatically  as  the  waypoints  are  overflown.  A 
software  check  is  implemented  such  that  the  pilot  must  rotate  the  OBS  knob  to 
the  appropriate  course  before  the  active  waypoint  number  is  incremented.  The 
bottom  rotary  switch  controls  the  initial  start-up  mode  of  the  GPS  receiver. 

In  Auto  Start  mode,  the  system  assumes  that  the  GPS  receiver  has  not  moved 
substantially  since  the  last  time  the  GPS  set  was  turned  on.  In  this  case, 
the  system  uses  the  position  fix  stored  in  non-volatile  memory  from  the  last 
flight  as  the  assumed  initial  location  of  the  receiver.  In  the  Manual  Start 
mode,  the  operator  oust  enter  the  initial  position  of  the  receiver  to  the 
nearest  degree  in  latitude  and  longitude. 

The  CDU  design  is  based  on  the  Z80  microprocessor,  and  employs  the  STD 
BUS  form  factor  and  bus  structure. 

3.8  Ins  t  rumentat ion 

The  GPS  T  and  E  equipment  was  equipped  with  instrumentation  to  record: 

a)  receiver  channel  and  position  software  activity  and  outputs, 

b)  received  carrier  to  noise  ratios  on  each  satellite, 

c)  pilot  (CDU,  CDI,  OBS)  activity,  and 

d)  aircraft  attitude  (roll,  pitch,  heading). 

Receiver  channel  and  position  software  activity  is  monitored  by  recording 
inter-task  messages.  Table  3-13  shows  the  general  content  of  the  various 
messages  recorded  on  magnetic  tape  during  a  typical  flight.  Tables  3-14 
through  3-21  show  the  message  formats  frequently  used  in  assessing  status  and 
position  estimation  performance.  Table  3-22  defines  the  alrcrsft  attitude 


TABLE  3-11 


ACTIVE  WAYPOINT  DATA  FILE 


NAME 

DESCRIPTION 

ACCURACY 

KEY 

UPDATE  INTERVAL 

PLAT 

Present  Latitude 

o.r 

POS 

2 

sec 

PLNG 

Present  Longitude 

0.1’ 

POS 

2 

sec 

GS 

Ground  Speed 

l.Okt 

GS 

2 

sec 

nt 

Actual  Track  (Magnetic) 

1.0  deg 

GS 

2 

sec 

WPT 

Active  Waypoint 

N/A 

WPT 

2 

sec 

DIST 

Distance  to  Active  Waypoint 

0.1  nm 

R/B 

2 

sec 

BRG 

Bearing  to  Active  Waypoint 

1.0  deg 

R/B 

2 

sec 

DTK 

Desired  Track 

1.0  deg 

CRS 

2 

sec 

NIK 

Next  Desired  Track 

1.0  deg 

CRS 

2 

sec 

TTG 

Time-to-go 

0.1  mLn 

TTG 

2 

sec 

XTK 

Cross-track  Deviation 

0.01  nsd. 

XTK 

2 

sec 

TKE 

Track  Angle  Error 

0.1  deg 

XTK 

2 

sec 

VAR 

Magnetic  Variation 

1.0  deg 

GMT 

2 

sec 

GMT 

Greenwich  Mean  Time 

1  sec 

GMT 

1 

sec 

OBS 

Current  OBS  Setting 

1  deg 

OBS 

1 

sec 

CSE 

Course  Selection  Error 

0.1  deg 

OBS 

1 

sec 

TDIST 

Total  distance  flown 
from  Waypoint  0. 

1.0  nml 

TTG 

2 

sec 

KEY 

DISPLAY  (Note  1) 

EXPLANATION 

R/B 

WPT 

RNG  19.2 

NMI 

Range  to  active  waypoint 

R/B 

WPT 

BRG  207 

MAG 

Bearing  to  active  waypoint 

TTG 

WPT 

TTG  12.1 

MIN 

Time  to  go  to  active  waypoint 

TTG 

TOT 

TIME  2:31 

HRS 

Total  time  in  air  from  waypoint  0. 

XTK 

XTK 

DEV  0.02 

NMI 

Cross-track  deviation 

XTK 

TRK 

ERR  14.1 

DEG 

Track  Angle  Error  (DTK-TK) 

CRS 

CRS 

TO  205 

MAG 

Course  to  active  waypoint  (Note  2) 

CRS 

CRS 

FROM  241 

MAG 

Course  from  active  waypoint  (Note  2) 

GS 

GS 

135 

KTS 

Ground  Speed 

GS 

TK 

206 

MAG 

Actual  Track 

GMT 

GMT 

1453:11 

HRS 

Greenwich  Mean  Time 

GMT 

MAG 

VAR  15 

DEG 

Magnetic  Variation 

POS 

LAT 

42DG 

17.0* 

Present  Latitude 

POS 

LONG  7 IDG 

43.6' 

Present  Longitude 

OBS 

OMNIBRG  170 

DEG 

OBS  Setting 

OBS 

OBS 

ERR  2.1 

DEG 

Course  Selection  Error  (DTK  -  OBS) 

VOR 

***INVALID  ENTRY 

Not  used  in  OPR  Mode 

INT 

***INVALID  ENTRY 

Not  used  in  OPR  Mode 

L/L 

WLAT  4 IDG  23 

.6' 

Active  waypoint  latitude 

L/L 

WLNG  72DG  42 

.1' 

Active  waypoint  longitude 

AWP 

ART 

WP:  GDM 

VOR 

Reference  ldent  and  type 

AWP 

Ill. 

.OMAG  16.0 

NMI 

Bearing  and  range  from  reference 

Note  1:  The  first  display  listed  for  each  key  appears  the  first  tine  the  key  is 

struck,  (except  CRS  key,  see  Note  2.)  Subsequent  keystrokes  cause  the  two 
display  lines  to  alternate. 

Note  2:  The  first  strike  of  CRS  Key  causes  "CRS  TO"  to  be  displayed  if  current 
direction  is  "TO"  active  waypoint.  Otherwise,  "CRS  PROM"  is  displayed. 
Subsequent  keystrokes  cause  lines  to  alternate. 


TABLE  3-13 


DATA  RECORDING  MESSAGES 


Message 

Block 

Number 

Source 

Content 

Frequency 

0 

DAU 

System  Reset 

Startup 

10 

DAU 

Receiver  Configuration  Data 

Startup 

15 

NSW 

PNP-DAU  Time  Synchronization 

Startup 

20 

PSW 

Position  Estimate  Sent  to  NSW 

l/sec 

25 

DAU 

Operator  Startup  Data  (if  defined) 

Startup 

30 

PSW 

Raw  GPS  Satellite  Nav  Data  (1  word) 

Variable 

40 

PSW 

Acquisition  Command 

As  Required 

50 

PSW 

Phased  Transition  Command 

As  Required 

60 

PSW 

Navigation  Command 

As  Required 

70 

PSW 

Measurement  Command  Data 

1-5/ sec 

80 

PSW 

Satellite  Reassignment  Command 

As  Required 

90 

PSW 

Nav  Plus  Data  Command 

Single  Channel 

110 

PSW 

Ephemeric,  Clock  and  Almanac 

Inf  requent 

120 

PSW 

PSW  Error  Message 

130 

NSW 

Receiver  Startup  Data  to  PSW 

Startup 

140 

PSW 

PNP  Status  to  DAU 

3/Mlnute 

150 

PSW 

Raw  Position  Fix 

1/2.2  sec  Typical 

160 

PSW 

Tracked  Position  Fix 

l/sec 

170 

NSW 

Pilot  Flight  Plan  Waypoint  Data 

As  Required 

180 

NSW 

CDU,  CDI,  OBS  Data 

l/sec 

190 

NSW 

Pilot  Startup  Data  From  CDU 

Startup 

200 

NSW 

Displayed  Waypoint  Data  Type 

When  Changed 

210 

NSW 

OBS  Data 

l/sec 

220 

NSW 

Displayed  Waypoint  Data 

l/sec 

230 

NSW 

Clock  Status  Sent  to  CDUI 

Startup 

400 

ADU 

RAS  Data 

l/sec 

410 

ADU 

DME  Data 

Variable 

420 

ADU 

ADU  Status 

Variable 

430 

ADU 

ADU  Error 

Variable 

DAU  Data  Acquisition  Unit 
NSW  Navigation  Software 
FSW  Position  Software 

ADU  Aircraft  Data  Unit 

RAS  Real  Tine  Aircraft  State 
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TABLE  3-14 


BLOCK  #10  RECEIVER  CHANNEL  CONFIGURATION  EATA 


Definition 

Format 

Units 

1 

Block  ID  (-10) 

constant 

N/A 

2 

Block  Length  (-14) 

constant 

words 

3 

Dual/Single  Channel 

binary,  unsigned 

N/A 

4 

5-6 

-  0  Dual 

-  377g  Single 

Statistics  Collect  (1  -  collect) 

AGC  Loop  Corner  Frequencys 

byte 

N/A 

7-8 

Nominal  AGC  Gains 

byte 

N/A 

9-12 

Search  Thresholds 

integer 

N/A 

13 

Re-search  Frequency  offset 

integer 

Hz 

14-16 

AFC  Loop  Noise  Band widths 

byte 

Hz 

17-18 

Code  Loop  Noise  Band widths 

byte 

Hz 

19-26 

AFC  Lock  Gains 

integer 

N/A 

27-28 

Bit  Sync  Histogram  Threshold 

byte 

N/A 

TABLE  3-15 

BLOCK  #30  RAW  GPS  NAV  MESSAGE  DATA  WORD 


Byte 

Definition 

Format 

Units 

1 

Block  ID  (-30) 

constant 

n/a 

2 

Block  Length  -  (13) 

constant 

bytes 

3-10 

GPS  Time 

real 

Becond  s 

11-14 

Mission  Time 

Integer 

mi  Hi  sec  s 

15 

Channel  Status 

see  block  #70 

16 

SV  ID 

binary,  unsigned 

N/A 

17 

Subframes  ID 

6  LSBs 

binary,  3  LSBs 

N/A 

18 

Word  count  (1-10)  within  subframe 

binary,  4  LSBs 

N/A 

19 

Parity  check  flag:  0  -  data  valid 

binary,  LSB 

N/A 

20 

1  -  parity  fail 

Least  significiant  byte  of  demodulat- 

bl nary 

N/A 

21 

ed  parity  stripped  data  word 

2nd  byte  of  demoudlated,  parity 
stripped  data  word 

binary 

N/A 

22 

Most  significant  byte  of  demodulated 
parity  stripped  data  word 

binary 

N/A 

23-24 

Wide  Band  Power 

see  block  #70 

25-26 

Narrow  Band  Power 

see  block  #70 

TABLE  3-16 

BLOCK  #40  ACQUISITION  COMMAND  DATA 


mi 

Defi nit ion 

Format 

Units 

LSB  Scale 

i 

Block  ID  (-40) 

constant 

N/A 

1 

2 

Block  Length  (-12) 

constant 

bytes 

1 

3-10 

GPS  Time 

real 

second  s 

- 

11-14 

Mission  Time 

binary,  unsigned 

milli¬ 
second  8 

20  ms. 

15 

SV  ID  Assignment  (1-31) 

binary,  unsigned 
6  LSB's 

N/A 

1 

16-17 

Carrier  frequency  prepositioning 
estimate  (Note  1) 

16  -  LSByte 

17  -  MSByte 

binary,  2's 
complement 

Hz 

4.064 

18 

Sub-chip  phase  code  prepositioning 
estimate  (0-9) 

binary,  unsigned 
4  LSB'S 

C/A 

code  chip 

0.1 

19 

Chip  phase  code  prepositioning 
estimate,  (0-1022)  least 
significant  byte 

binary,  unsigned 

C/A  code 
chip 

1 

20 

Chip  phase  code  prepositioning 
estimate,  2  MSB's 

binary,  unsigned, 
2  LSB's 

C/A  code 
chip 

256 

21 

Epoch  phase  code  prepositioning 
estimate  (0-19) 

binary .unsigned , 
6  LSB's 

C/A 

epochs 

1 

22 

Search  and  Acquisition  Strategy 

(Note  2) 

(Note  2) 

(Note  2) 

23-24 

Bit  Count  (0-299)  (Note  3) 

binary,  unsigned 

quanta 

1 

23  -  LSByte 

24  -  MSByte 

binary,  unsigned 

1  LSB 

quanta 

256 

TABLE  3-16.  (CONT'D) 


(1)  All  prepositioning  estimates  are  valid  at  the  trailing  edge  of  the 
particular  20  msec  system  clock  pulse  defining  the  beginning  of 
acquisition  of  an  SV. 

(2)  Search  and  Acquisition  Strategy  Byte,  contains: 

-  Bit  7  (MSB)ls  DATA  DEHOD  FLAG 

0:  channel  to  demodulate  data 
1:  channel  not  to  demodulate  data 

-  Bit  6  is  SYNC  FLAG 

0:  channel  to  perform  Bit/Subframe  Sync  operations 
i:  channel  to  use  PNP  given  BIT  COUNT  (Bytes  10  and  11) 
to  drive  Subframe  Sync 

(DATA  DEMOD  FLAG  -  1  SYNC  FLAG  -  0  is  Illegal  combination) 

-  Bit  5-0  is  HALF  SEARCH  APERTURE  CODE  (HSAC) 

(format:  binary,  unsigned)  related  to  half  search  aperture  as 
follows 

HSAC 
~0 

1  16 

2  24 


63  512 

Half  Search  Aperture  (chips) 

l.e. ,  HSAC  -  - 1 

8 

(3)  Bit  Count  is  valid  only  if  SYNC  FLAG  ■  1 


Ts7 


. 


TABLE  3-17 

BLOCK  #70  MEASUREMENT  DATA 


Byte 


1 

2 

3-10 

11-14 

15 

16 

17 

18 

19 

20 
21 
22 
23 

24-25 

26-27 

28 

29 


Definition 

Format 

LSB  Scale 

mmmM 

Block  ID  (-70) 

constant 

N/A 

1 

Block  Length  (-15) 

constant 

words 

1 

GPS  lime 

real 

seconds 

— 

Mission  Time 

Channel  Status  (Note  2) 

binary,  unsigned 

milli¬ 
second  B 

20  ms. 

SV  ID  (1-31) 

binary,  unsigned, 
6  LSB's 

N/A 

1 

Least  significant  byte  of  receiver 
channel  estimated  carrier  frequency 

binary  2'coraple- 

ment 

mHz 

6.35 

(Note  2) 

Second  byte  of  receiver 

channel  estimated  carrier  frequency 
Most  significant  byte  of  receiver 
channel  estimated  carrier  frequency 

Sub-chip  code  phase,  valid  at  dwell 

binary,  unsigned, 
4  LSB's 

C/A 

0.1 

ending  (0-31) 

code  chip 8 

Chip  code  phase,  least  significant 
byte.  Valid  at  dwell  ending  (0-1022) 
Chip  code  phase,  2  MSB's 

binary,  unsigned 

C/A 

code  chips 

1 

binary,  unsigned, 
2  LSB's 

C/A 

code  chips 

256 

Epoch  code  phase  valid  at  dwell 

binary,  2'comple- 

C/A 

1 

ending  (-127  to  +128)  (Note  3) 

ment 

code  epoch 

Wideband  Power 

24  -  LSByte 

binary,  unsigned 

quanta 

25  -  MSByte 

Narrow  Band  Power 

binary,  unsigned 

quanta 

26  -  LSByte 

27  -  MSByte 

N/A 

Data  word  count  (0-9)  within  sub- 
frame  at  dwell  ending 

bl nary ,  unsigned , 
4  LSB 

1 

Data  bit  count  (0-29)  within  sub- 
frame  at  dwell  ending 

binary,  unsigned, 
5  LSB 

N/A 

1 

TABLE  3-17.  (CONT'D) 


Notes: 

(1)  During  navigation  node,  this  data  is  valid  at  the  leading  edge  of 
the  particular  20  nsec  system  clock  pulse  that  defines  satellite 
transition. 


(2)  Channel  status  byte 


|7  6  5  4  |  3  2  1  f0| 


♦Bit  0 
Bit  1 

♦♦Bit  2 
Bit  3 
Bits  4-7 


Set  when  receiver  fault  sensed 

Init  Node  or  Configuration  Node  Complete,  Carrier  Lock, 
or  Neasurement  Valid 
Frame  Synch  Complete/ In  Data  Demod 
1  Pass  of  Search  Completed 
Receiver  State  (coded) 


0  -  Idle 

1  *  In  Limited  Search  (±  15  chip) 

2  ■  In  i  2  chip  Search 

3  ■  In  ±  1  chip  Search 

4  *»  In  False  Alarm  Check 

5  ■  In  Re-search 

6  ■  In  AFC  Pull-in  #1  (Pre  Bit  Sync) 

7  -  In  AFC  Pull-in  #2  (Pre  Bit  Sync) 

8  -  AFC  Lock 

9  «  In  AFC  Pull-in  #1  (Post  Bit  Sync) 

10  -  In  AFC  Pull-in  #2  (Post  Bit  Sync) 

11  *  Fine-Time 

12  -  Bit  Synch  Operations 

13  -  Frame  Synch  Operations 

14  -  Data  Demodulation 

15  “  Statistics  Done 


(3)  Pseudo-range  can  roll  over  a  bit  edge  (overflow  or  underflow), 
therefore  it  can  be  larger  than  20msec,  or  it  can  be  negative; 
referenced  to  the  initial  position  value. 


(4)  Hide  band  (WBP)  and  narrowband  (NBP)  power  quantities  are  used  in 
the  following  equation  to  compute  the  C/N0  estimate: 

(NBP  -  WBP) 

C/N0  -  10  log10  2000  -  dB  Hz 

(40  WBP  -  NBP) 


Note:  6  dB  Mist  be  added  to  C/N0  estimate  if  data  is  not  being 
demodulated. 


♦Pseudo-Range  is  valid  if  Bit  1  is  set  in  Acquisition,  Transition  or 
Navigation  Phase.  Bit  1  is  set  when  mode  is  complete  when  in 
Initialisation  or  Configuration  Nodes. 

♦♦Bit  2  is  set  when  frame  synch  is  achieved  and  data  is  being 
demodulated,  unless  Frame  Sync  is  not  required.  Then  it  indicates 
data  demodulation. 
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WORD 


DEFINITION 


FORMAT 


UNITS 


LSB  SCALE 


1 

Block  ID  («  200) 

constant 

- 

1 

2 

Block  Length  («14) 

constant 

words 

1 

3-4 

Mission  Time 

positive  integer 

milliseconds 

20  ms. 

5(bits  8-15) 

Active  Waypoint  number 
(0  to  19) 

integer 

- 

1 

5(bits  2,1,0) 

CDI  Sensitivity 

integer 

- 

1 

0  “  Enroute 

1  -  Terminal 

2  -  Approach 

6-9 

Active  Waypoint  Latitude 

double  real 

radians 

- 

10-13 

Active  Waypoint  Longitude 

double  real 

radians 

- 

14(bits  0-5) 

RNAV  Parameter  Number 
(11  through  32) 

11  -  Range-to-Waypoint 

12  -  Bearing-to-Waypoint 

13  -  Time-to-go 

14  -  Not  used 

15  ■  Cross-track  deviation 

16  -  Track  angle  error 

17  «  Course  to  waypoint 

18  m  Course  from  waypoint 

19  ■  Ground speed 

20  -  Track  angle 

21  -  Greenwich  Mean-Time 

22  ”  Not  used 

23  -  Present  latitude 

24  ■  Present  Longitude 

25  ■  OBS  setting 

26  •  Not  used 

27  ■  Artificial  Waypoint 

reference  type 

28  ■  Artificial  waypoint 

reference  ident. 

29  ■  Artificial  waypoint 

offset  range 

30  *'  Artificial  waypoint 

offset  bearing 

31  «  Waypoint  latitude 

32  «  Waypoint  longitude 

integer 

1 

TABLE  3-21 


i 


AREA  NAVIGATION  DATA  BLOCK  #220  (TO  CDU) 


WORD 

DEFINITION 

FORMAT 

UNITS 

LSB  SCALE 

1 

Block  ID  (-220) 

constant 

i 

i 

1 

2 

Block  length  (-20) 

constant 

16-bit  words  | 

1 

3-4 

Mission  Time 

positive  integer 

Billiseconds 

20  as. 

5(blt  8) 

NAV  flag  0  -  okay 

1  -  not  OK 

Integer 

- 

5(bits  1,0) 

To/From  flag  0  *  from 

1  -  to 

2  -  neither 

integer 

— 

• 

6(bitsl0,9,8) 

Solution  Quality  0  -  OPR 

1  -  UNR 

2  -  LCH 

3  -  BAD 

4  -  STR 

integer 

6  (bit*  0*5) 

RNAV  Paraaeter  No.  (see  block 
#200) 

integer 

- 

1 

7-8 

RNAV  Paraaeter 

real 

— 

9-10 

OBS  Setting 

real 

radians 

- 

11-12 

CDI  Deflection 

real 

dots 

1 

13-16 

Present  Latitude 

double  real 

radians 

- 

17-20 

Present  Longitude 

double  real 

! 

radians 

1 

8 


TABLE  3-22 


REAL-TIME  AIRCRAFT  STATE  DATA 


Parameter 

Value  of 

LSB 

Range 

Min 

Max 

Heading 

(360V512) 

0 

360°  (binary  512) 

Roll 

1° 

-31°  (binary  0) 

+32°  (binary  63) 

Pitch 

1° 

-31 #  (binary  63) 

+32°  (binary  0) 

Outside  air 
Temperature 

0.625°C 

-40#C  (binary  0 

+39.375°C  (binary  127) 

Indicated 

Air  Speed 

1  mm  H20 

0  (binary  0) 

1023  (binary  1023) 

Rate  of 

Climb 

24  mm  ^O/Mn 

-768  (binary  0) 

+744  (binary  63) 

Mote:  Pitch  reading  is  not  linear  for  pitch  angles  greater  than  20  degrees. 
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The  receiver  carrier-to-noise-density  ratio  (C/N0)  ie  estimated  using  the 
wide-hand  energy  detected  by  the  AGC  loop  and  the  narrow-band  detector  shown 
in  Fig.  3-14.  Scaled  estimates  are  recorded  for  each  measurement  in  blocks  30 
and  70.  C/N0  is  then  calculated  using  the  formula  in  Note  (4),  Table  3-17. 

Aircraft  altitude  is  measured  using  standard  aircraft  Instruments  which 
have  been  connected  via  a  digital  interface  to  a  Z80-based  Aircraft  Data  Unit 
assembly  (ADU).  The  ADU  in  turn  communicates  via  a  serial  RS232  interface  to 
the  DAU,  nominally  once  per  second.  The  altitude  parameters  are  defined  in 
Table  3-19. 

Pilot  activity  is  monitored  by  recording  request  and  response  messages 
between  the  DAU  and  CDU  using  the  formats  shown  in  Tables  3-20  and  3-21. 


4.0  MEASURED  PERFORMANCE 


4.1  Test  Facilities 


4.1.1  Ground  Test  Laboratory 


During  Integration  of  the  GPS  T  and  E  system  a  tower-mounted  antenna  was 
employed.  Initial  evaluation  of  the  hardware  and  software  functions  were 
completed  in  the  laboratory  before  flight  tests  commenced. 


4.1.2  Aircraft  Installtlon 

The  GPS  T  and  E  equipment  was  Installed  in  a  Rockwell  Aerocommander  300  A, 
as  shown  in  Figs.  4-1  through  4-4.  The  T  and  E  equipment  was  configured  as 
shown  in  Figs.  4-5  and  4-6.  The  CDU  and  CDI  are  installed  as  shown  in  Fig. 
4-7.  The  traditional  difficulty  of  pushing  buttons  on  a  vertically  mounted 
cockpit  device  was  relieved  for  the  CDU  because  the  engine  controls  provide  a 
convenient  stabilizing  hand  rest  for  most  power  settings.  The  performance 
envelope  for  the  Aerocommander  is  shown  in  Table  4-1. 


4.1.3  Mode  S  Experimental  Facility 


Position  truth  was  established  using  the  Mode  S  Experimental  Facility 
(MODSEF)  operated  at  Lincoln  Laboratory.  It  includes  a  monopulse  surveillance 
system  capable  of  Interrogating  ATCRBS  and  Mode  S  Transponders.  Its  principal 
characteristics  are  provided  in  Table  4-2. 


During  GPS  tests  Mode  S  Experimental  Facility  surveillance  accuracy  was 
verified  using  a  calibration  transponder  located  6  miles  from  the  Mode  S 
Experimental  Facility.  The  aircrafts  altimeter  calibration  was  also  verified. 
Since  most  of  the  Hanscom  runways  are  visible  to  the  Mode  S  Experimental 
Facility  sensor,  a  test  was  conducted  in  which  the  aircraft  parked  at  a  taxi 
way  reference  point,  to  verify  surveillance  range/altitude. 

The  Mode  S  Experimental  Facility  surveillance  data  for  the  GPS  test 
aircraft  was  recorded  and  identified  based  on  either  an  ATCRBS  discrete  code 
or  Mode  S  ID  (both  transponders  were  available  on  the  Aerocommander). 
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-  PA8SENGER 
COMPARTMENT 


M0010  001 


7ZC 


P130-1683 


■racaa 


Complete  System 

Without  Instrumentation 

Cruise  Speed 

160  KIAS 

160  KIAS 

Celling 

10  Kft. 

10  Kft. 

Hue  at  10  Kft. 

1.5  Hours 

2.5  Hours 

Assumptions 

2  Pilots 
1  Operator 


Antenna 

4*  monopulse 

Transmit 

150  watts  with  SLS 

Receiver 

Phase  monopulse 

10  dB  noise  figure 

RSLS 

Nominal  Range 

75  nmlles 

Surveillance  Accuracy 

Range 

40  feet  lo 

Azimuth 

0.05  deg  lo 

Calibration 

Test  transponder  at  surveyed  site 

6  miles  north  of  The  Mode  S 
Experimental  Facility 

Display 

TPX  42  display 

4.1.4  Analysis  Software 


The  data  recorded  on  the  aircraft,  and  at  the  Mode  S  Experimental 
Facility  were  processed  as  shown  in  Fig.  4-8  to  produce  a  basic  surveillance 
data  base  and  error  statistics. 

The  difference  between  the  GPS  and  Mode  S  Experimental  Facility  position 
estimates  was  determined  by  using  every  ungarbled  Mode  S  Experimental  Facility 
report  (unsmoothed,  occurring  at  4  second  rate),  and  the  two  adjacent  GPS 
position  estimates  as  shown  in  Fig.  4-9.  Conversion  from  the  Mode  S 
Experimental  Facility  R,  8,  H  position  estimate  to  the  aircraft's  XYZ  in  earth 
centered  earth  fixed  ECEF  coordinates  was  done  using  matrix  transformations* 
which  account  for  the  WGS-72  ellipsoid  of  revolution  earth  model. 

Plots  of  XT  versus  time,  GPS  position  error  with  respect  to  the  Mode  S 
Experimental  Facility  versus  time,  and  error  statistics  were  developed. 

GDOP  was  independently  estimated  using  Z  STARS,  a  Lincoln  program  which 
provides  visibility  to  current  satellites.  Using  Z  STARS,  GDOP  was  calculated 
based  on  all  satellites  actually  in  track. 

4.2  Laboratory  Tests 
4.2.1.  Antenna 

The  Chu  Associates,  Inc.,  model  CA-3224  volute  antenna  was  tested  while 
mounted  on  a  5-foot  by  5-foot  aluminum  plate  shaped  to  simulate  the  local  Aero 
Commander  fuselage  as  shown  in  Fig.  4-10.  Anecholc  chamber  azimuth  gain 
measurements  were  made  at  various  elevation  angles.  The  results,  shown  in 
Figs.  4-11,  -12  and  -13,  show  that  the  specification  was  met  or  exceeded  over 
most  elevation  angles  except  within  the  range  5-10*  where  the  gain  was  2-dB 
lower  than  specified.  Measurements  above  50*  elevation  were  not  possible  due 
to  mechanical  conflicts  between  the  fuselage  plate  and  the  anechoic  chamber 
mount.  It  is  assumed  that  a  gain  of  0  ±  1  dBl  in  the  region  50*-90*  has  been 
provided. 

The  elevation  gain  profile.  Fig.  4-14,  was  developed  by  computing  the 
mean  azimuth  gain  at  each  elevation  angle  from  samples  taken  every  10*  in 
azimuth.  The  standard  deviations  in  azimuth  gain  were  also  computed,  as  shown 
by  the  dotted  lines. 


*  See  Appendix  C  for  details. 
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4.2.2  Receiver  Channel  Performance 


The  receiver  factory  acceptance  tests,  conducted  using  the  equipment 
configuration  shown  in  Fig.  4-15,  verified  compliance  with  all  receiver 
specifications  including: 

a)  Receiver  Hardware  -  PNP  Software  Interface. 

b)  Acquisition  tine  versus  doppler  and  C/Hq. 

c)  Data  error  rate  versus  C/H0> 

d)  Transition  performance  versus  C/N0. 

e)  Navigation  perforaance  versus  C/N0. 

f)  Nav  plus  Data  perforaance  versus  C/N0. 

g)  AFC  error  perforaance. 

h)  Pseudo-range  accuracy. 

i)  C/N0  aeasureaent  accuracy. 

During  the  acceptance  tests,  acquisition  of  all  satellite  codes  now  in 
use  with  doppler  offsets  up  to  >  1000  Hz  (>  750  Hz  requireaent)  was 
demonstrated.  Tests  of  the  Transition,  Navigation  and  Navigation  plus  Data 
modes  were  accomplished  by  phase-locking  the  receiver  oscillator  to  the 
satellite  simulator  oscillator  and  manually  entering  prepositioning  code  and 
doppler  data  into  the  PNP  simulator.  CPU  activity  aonltorlng  verified  that 
the  receiver  Z8000  microcomputer  task  execution  times  did  not  exceed  their 
allotted  times  in  any  mode.  Finally,  all  visible  satellites  were  successfully 
acquired  using  a  roof  antenna.  All  tests  were  conducted  on  both  channels. 

Following  shipment  of  the  receiver  to  Lincoln  Laboratory,  the  acceptance 
tests  were  repeated  in  a  laboratory  using  the  same  test  configuration  as  in 
Fig.  4-15  except  that  a  volute  antenna  identical  to  that  installed  on  the 
Aerocommander  was  used. 

Table  4-3  summarizes  the  major  requirements  and  performance  measured 
during  the  factory  acceptance  tests  and  during  the  post-acceptance  evaluation. 
Results  of  specific  measurements  are  shown  in  Figs.  4-16  through  4-24,  and  the 
results  of  a  test  in  which  a  GPS  satellite  was  tracked  for  5  hours  are  shown 
in  Fig.  4-25. 

The  measured  receiver  performance  was  also  evaluated  with  respect  to 
theoretical  predictions.  The  results,  shown  in  Figs.  4-26,  -27,  and  in 
Table  4-4,  show  good  agreement. 

To  assess  the  5-hour  tracking  test  data,  the  expected  receive  C/Nq  was 
calculated  as  shown  in  Table  4-5.  Worst-case  RF  power  (C/A  code  only)  at  the 
user  antenna  is  seen  to  be  -160  dBw.  Note  that  the  receiver  preamp  is  assumed 
to  be  located  at  the  antenna  for  minimum  cabling  loss. 
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FIG.  4-15.  DUAL  CHANNEL  RECEIVER  ACCEPTANCE  TEST  CONFIGURATION 


Parameter 


Requirement 


Measured  Performance 


Acquisition 

Time  to  acquire 

AFC  errors 

Bit  Error  Rate 

C/Nq  Error 

25  sec 

(based  on  6  min  TTFF) 

Aj  <  10  Hz  for  data 
demodulation 

<  10“5 

±  1  dB 

25  sec,  P  >.9 

of  <  1.4  Hz  for  C/Nq 

>  36  dBHz 

<  10-5 

°C/No  <  °*7  dB 

Phased  Transition 

AFC  Error 

Not  specified 

af  <  5  Hz 

Pseudo  Range  Error 

Not  specified 

0T  <  45  ns  for  C/NQ 

>  36  dBHz 

C/N„  Error 

±  1  dB 

°C/No  *  *85  dB 

Navigation 

AFC  Error 

Not  specified 

Oj  <  7  Hz 

Pseudo  Range  Error 

oT  <  50  ns 

oT  <  47  ns  for  C/NQ 

>  36  dBHz 

C/N0 

±  1  dB 

aC/No  *  1,0  dB 

Nav  Plus  Data 

Same  as  Navigation 

Same  as  Navigation 

Same  as  Navigation 

and  Acquisition 

and  Acquisition 

and  Acquisition 

Acquisition 

4.5  dB  at  5*  elevation. 

Link  Margin* 

At  least  4  dB 

10  dB  at  40*  elevation, 
for  the  volute  antenna. 

RECEIVER  INPUT  C/No  dBHz 
RECEIVER  AFC  PERFORMANCE  IN  ACQUISITION  PHASE 


RECEIVER  INPUT  C/No  dBHz 

RECEIVER  C/No  ESTIMATOR  PERFORMANCE  IN  ACQUISITION  PHASE  r42-2l90 


RECEIVER  INPUT  C/No  dBHz 

FIG.  4-22.  RECEIVER  PSEUDO  RANGE  PERFORMANCE  IN  THE  NAVIGATION  STATE 


RECEIVER  IN  C/No  dBHz 

RECEIVER  AFC  PERFORMANCE  IN  THE  NAVIGATION  STATE 


CH  1 


RECEIVER  INPUT  C/No  dBHz 

FIG.  4-24.  RECEIVER  C/No  ESTIMATOR  PERFORMANCE  IN  NAVIGATION  STATE 


RECEIVER  PERFORMANCE  DURING  SATELLITE  PASS 


STATE 


FIG.  4-26.  COMPARISON  OF  THEORETICAL  AND  MEA8URED  RECEIVER  FREQUENCY  TRACKING  PERFORMANCE 


TABLE  4-4 


COMPARISON  OF  THEORETICAL  AND  MEASURED  RECEIVER  PERFORMANCE 


Characteristic 

JggggJj 

Minimum  Acquisition  C/N0 

BIF  -  1000  Hz 

BL  -  10  Hz 

L  -  2  dB 

35.7  dB  Hz 

35  dB  Hz 

Nav  Loss  of  Lock 

BIF  -  200  Hz 

BL  -  10  Hz 

L  -  2  dB 

31.5  dB  Hz 

33  dB  Hz 

TABLE  4-5 

USER  RECEIVED  POWER  FOR  C/A  CODE  ONLY 


Satellite  Transmitter 
Power  (dBw) 

RF  Loss  (dB) 

Antenna  Polarisation 
Loss  (dB) 

Antenna  Gain 


Satellite 
Overhead  (at 
Zenith) 


14.25 

1.0 

0.25 

15.0 


Satellite 
at  Elevation 
Angle  of  5* 

14.25 

1.0 

0.25 

12.0 


Satellite  ERP  (dBw) 

Path  Loss  (dB) 

Atmospheric  Absorption  (dB) 

28.0 

182.5 

25.0 

184.2 

0.8 

Power  at 

User  Antenna  (dBw) 

-154.5 

-160.0 

Receiver  noise  (d Bw 

(4  dB  noise  Figure) 

-200 

-200 

Receiver  C/No  (dBHz) 


45.5 


40.0 


It  Is  evident  from  Fig.  4-25  that  the  measured  C/N0  values  were  greater 
than  expected.  The  mean  C/N0  for  elevational  angles  greater  than  15*  was 
47.8  dB-Hz,  2.3  dB  greater  than  the  45.5  dB-Hz  prediction  for  a  zenith  aspect 
It  was  postulated  that  the  additional  power  was  due  to  excess  satellite 
transmit  power,  excess  receiver  antenna  gain,  and  excess  sensitivity  in  the 
receiver  preamplifier. 

To  verify  this  a  calculation  of  the  expected  received  C/N0  and  link 
margin  referred  to  the  preamplifier  input  was  made,  as  shown  in  Table  4-6,  and 
plotted  in  Fig.  4-25.  It  is  apparent  that  for  satellite  number  6  the  link  has 
about  2  dB  additional  margin,  due  probably  to  excess  satellite  radiated  power. 
Also,  the  variations  in  C/N0  were  examined  and  found  to  be  generally 
consistent  with  the  ripple  in  the  volute  antenna  patterns. 

4.2.3  Static  Acquisition  Performance 

Initial  laboratory  tests  using  the  tower  mounted  volute  antenna  indicated 
that  the  T  and  E  Equipment  meets  the  Time  to  First  Fix  Requirements  of 
6  minutes  maximum  when  the  receiver  is  initialized  with  current  almanac  data. 

The  sequence  of  events  in  a  representative  static  acquisition  is  shown  in 
Table  4-7. 

4.2.4  Static  Position  Accuracy  and  Stability 

Tests  were  conducted  to  characterize  the  static  accuracy  and  stability 
performance.  Test  data,  obtained  using  the  tower-mounted  antenna,  was 
collected  and  analyzed.  Results  from  a  representative  26-minute  test 
conducted  on  19  November  1982  using  five  GPS  satellites  were: 

a.  RMS  horizontal  error:  93  feet 

b.  RMS  vertical  error:  104  feet 

A  scatter-plot  of  the  static  test  results  is  shown  in  Fig.  4-28. 

A  static  test  was  also  conducted  on  9  May  1983  using  a  later  version  of 
the  GPS  position  software.  This  version  of  the  position  software 
(Version  3.0)  included  the  capability  of  navigating  with  three  GPS  satellites 
and  the  baro-altimeter .  The  RMS  horizontal  error  for  this  run,  lasting 
63  minutes,  was  95  feet,  nearly  the  same  as  the  five-satellite  case.  Although 
this  version  of  the  software  was  not  flight  tested,  the  basic  capability  for 
three-satellite  navigation  was  demonstrated. 


Mini  hub  user 
Receive  Power * 
dBw 


-160.0 

-159.7 

-159.0 

-158.0 

-158.8 

-160.0 


+1.50 

+1.74 

+1.86 

+1.92 

+1.94 

+1.95 


Correction 
for  Measured 
EDCR  Volute 
Antenna  Galn^ 
(dB) 

Net  C/N0 
(3  dB  Noise 
Figure)4 
(db  Hz) 

Link  Margin 
for  35  dBz 
Acquisition 
Threshold 
(db) 

-2.0 

40.5 

5.5 

0 

43.0 

8.0 

+1.9 

45.8 

10.8 

+1.2 

46.0 

11.0 

+0.4  (est) 

44.5 

9.5 

-0.4  (est) 

42.5 

7.5 

From  Rockwell  ICD-GPS-200,  May  81.  (Assumes  2  dB  atmospheric  loss). 

Digital  Communications  by  Satellite,  Spilker  Page  171. 

GPS-QTL-4-8 ,  page  10. 

The  nominal  noise  figure  for  the  EDCR  Is  4  dB.  The  actual  noise  figure  present 
during  the  characterisation  measurement  was  3  dB. 


TABLE  4-7 


STATIC  ACQUISITION  EXAMPLE 


Event 

Time 

(Minutes  :  Seconds) 

Start 

0 

Satellite  4  acquired 

0:32 

Satellite  6  acquired 

1:55 

Satellite  5  acquired 

2:49 

Satellite  8  acquired 

3:56 

First  Internal  Fix 

Transition  Mode  Started 

5:00 

Fix  Sent  to  Navigation  Software 

5:05 

Navigation  Mode  Started 

5:11 
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4.2.5  Interference  Effects 


Stetlc  tests  conducted  in  the  Aerocommander  Indies te  that  the  UHF  radios. 
VOR  receivers.  ATC  transponder  and  ME  did  not  appear  to  interfere  with  the 
operation  of  the  preaap  and  receiver  channels.  There  was  no  observed  effect 
on  the  operation  of  the  aircraft  avionics  when  the  GPS  receiver  was 
operated. 

In  addition,  a  s inula t ion  of  Interference  effects  in  the  Rockwell 
Aerocoaaander  test  aircraft  was  carried  out  for  the  FAA  by  the  DoD 
Electromagnetic  Compatibility  Analysis  Center  (ECAC).  The  ECAC  study  showed 
that  the  Aerocommander  avionics  would  not  interfere  with  the  operation  of  the 
GPS  receiver  for  the  particular  antenna  configuration  used  in  the  tests 
[Ref.  12]. 


4.2.6  Performance  Monitoring 

The  performance  monitoring  features  described  in  Section  3.5.6  were 
verified  during  static  tests.  Features  verified  Include  the  receiver  fault 
monitoring  via  receiver  flrmwave  diagnostics  and  receiver-position  processor 
interface  fault  detection  via  the  Position  Software. 

4.3  Flight  Tests 

During  the  course  of  the  GPS  test  and  evaluation  project,  a  total  of  25 
flights  were  made  during  which  a  total  of  33  hours  of  flight  data  was 
collected.  A  variety  of  flights  were  conducted  at  a  total  of  seven  different 
airports  as  listed  in  Table  4-8.  These  tests  measured  system  performance 
during  takeoffs,  landings,  approaches  and  turns.  The  tests  conducted  for 
airports  in  Massachusetts  and  southern  New  Hampshire  were  within  the  range  of 
MODSEF  surveillance. 

The  flight  test  program  was  divided  into  engineering  and  operational 
phases.  The  engineering  tests  were  conducted  primarily  to  determine  the 
effect  of  aircraft  dynamics  and  link  margin  on  the  accuracy  and  reliability  of 
the  GPS  navigation  system.  The  operational  tests  were  conducted  to  determine 
the  performance  of  the  system  under  a  variety  of  conditions  typical  of  general 
aviation  use  and  to  verify  the  compatibility  of  the  system  with  conventional 
air  navigation  systems. 

4.3.1  Engineering  Flight  Tests 

A  total  of  14  engineering  test  flights  were  made  during  which  a  total  of 
16.5  hours  of  data  was  collected.  These  flights  were  made  over  the  six-month 
period  lasting  from  June  to  December  of  1982.  Position  software  development 
continued  over  this  period,  so  that  most  of  the  effort  was  directed  at 
assessing  the  accuracy  of  the  system  in  various  dynamic  situations,  such  as 
level  flight,  turns,  climbing  and  descent.  At  the  beginning  of  December, 

1982,  the  position  software  was  frosen  at  version  2.1,  which  was  used  for  all 
subsequent  engineering  and  operational  flight  tests. 


TABLE  4-8 

AIRPORTS 

FOR  GPS  FLIGHT  TESTS 

Teat 

MODSEF 

Location 

Takeoff  Landing  Low-Approach 

Surveillance 

Hanacom  Field 

Lexington,  1 

MA 

Gardner  Municipal 

Gardner,  MA 

Logan  International 

Boston,  MA 

Bolre  Field 

Nashua,  NH 

Manchester  Airport 

Manchester, 

NH 

Burlington 

International 

Burlington, 

VT 

Dulles 

International 

Washington, 

DC 

4.3. 1.1  Engineering  Tests  with  Preliminary  Position  Software 

An  engineering  test  using  the  preliminary  position  software  was  conducted 
on  28  October  82.  The  results  of  this  test  are  included  here  to  illustrate 
the  benefit  of  tracking  all  satellites  in  view  and  to  quantify  the  fade  margin 
during  steep  turns. 

In  this  test,  the  aircraft  was  initially  in  a  holding  pattern  at  the 
LOBBY  intersection  as  shown  in  Fig.  4-29.  After  emerging  from  the  holding 
pattern,  the  aircraft  performed  two  360a  turns  at  a  bank  angle  of  30°.  The 
aircraft  next  was  flown  to  intercept  the  Shaker  Rills  waypoint  (SKR  NDB)  and 
then  landed  on  runway  11. 

Figure  4-30  shows  the  satellite  visibilities  during  the  28  October  1982 
flight.  All  five  operational  satellites  were  at  least  10°  above  the  horizon 
and  SV5  had  the  lowest  elevation  angle,  ranging  from  18°  to  25s  during  the 
run.  Figure  4-31  shows  the  GDOP,  HDOP  and  VDOP  values  during  the  flight.  The 
HDOP  and  VDOP  values  are  calculated  relative  to  the  local  MODSEF  coordinate 
system.  The  GDOP  value  ranged  from  S.l  to  5.4. 

During  the  flight,  the  aircraft  was  under  continuous  MODSEF  surveillance, 
allowing  the  system  navigation  accuracy  statistics  to  be  computed  by 
post-flight  processing.  The  position  accuracy  statistics  are  calculated  in 
two  ways  based  on:  1)  the  position  fixes  generated  every  2.2  seconds,  and  2) 
the  a-£  tracker  estimates  generated  once  per  second.  The  position  fix 
statistics  give  the  accuracy  of  the  Independent  position  fixes  in  the  normal 
navigation  mode.  The  tracker  estimate  statistics  give  the  system  accuracy  as 
displayed  to  the  pilot  using  linear  extrapolation  between  position  fixes. 

The  position  accuracy  statistics  for  the  28  October  flight  are  shown  in 
Table  4-9.  For  this  perllmlnary  version  of  the  position  software,  the 
horizontal  error  from  the  a-B  tracker  position  estimate  was  364  feet  (95Z). 
This  horizontal  error  value  easily  meets  the  AC90-45A  requirement  for 
non-precision  approach  and  is  better  than  the  500  ft  (95Z)  goal  for  straight 
and  level  flight.  The  rms  horizontal  error  was  211  feet,  which  also  meets  the 
FRP  requirements  for  non-precision  approach. 

The  portion  of  the  28  October  1982  flight  containing  the  30°  bank  angle 
turns  are  shown  in  Figs.  4-32a  and  4-32b.  As  shown  in  Table  4-10  the  95Z 
horizontal  position  estimate  error  was  429  feet,  a  factor  of  two  better  than 
the  1000  ft  95Z  accuracy  goal  during  30*  bank  turns.  Figures  4-33a  and  4-33b 
show  the  horizontal  and  vertical  errors  versus  time  during  the  30*  bank  turns 
for  the  position  fixes  and  the  trackes  estimates.  Figure  4-34  shows  the 
measured  carrler-to-nolse-density  (C/N0)  ratios  during  the  30*  bank  turns  for 
the  five  visible  satellites.  Note  that  SV5  drops  below  the  33  dB-Hz 
loss— of —lock  threshold  at  390  sec  and  515  sec.  These  two  occasions  correspond 
to  the  times  when  the  aircraft  was  banked  away  from  the  satellite.  From 
Fig.  4-30,  the  satellite  elevation  was  about  20*  so  that  the  satellite  was  at 
-10*  elevation  angle  with  respect  to  the  GPS  antenna. 
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TABLE  4-9 


POSITION  ACCURACY  STATISTICS  FOR  RUN  la.  28  OCTOBER  1982 


Position  Fixes: 


Error 

Mean 

Horizontal 

130 

Vertical 

185 

Tracker  Estimates: 

Error 

Mean 

Horizontal 

128 

Vertical 

194 

Standard 

Deviation 

RMS 

95Z 

111 

171 

276 

100 

210 

375 

Standard 

Deviation 

RMS 

95Z 

16'J 

211 

364 

126 

231 

423 

TABLE  4-10 


POSITION  ACCURACY  STATISTICS  FOR  30s  BANK  ANGLE  TURNS 


Position  Fixes: 

Error 

Mean 

Horizontal 

181 

Vertical 

166 

Tracker  Estimates : 

Error 

Mean 

Horizontal 

182 

Vertical 

191 

Standard 


Deviation 

RMS 

951 

54 

189 

245 

ft 

104 

196 

322 

ft 

Standard 


Deviation 

RMS 

95 X 

201 

271 

428 

ft 

135 

234 

423 

ft 

FIG.  4-33b.  POSITION  ESTIMATE  ERROR  DURING  30°  TURNS 


According  to  the  GPS  antenna  gain  characteristic  of  Fig.  4-14.  changing  the 
elevation  angle  from  +20*  to  -10*  causes  an  8  dB  drop  in  antenna  gain.  This 
loss  in  antenna  gain  would  drop  the  satellite  C/No  froa  the  aean  value  of  44.7 
dB-Hs  to  36.7  dB-Hz.  Since  the  C/No  value  for  SV5  dropped  to  less  than 
33  dB-Hz  during  portions  of  the  turn,  it  seems  likely  that  shielding  froa  the 
wings,  engine  nacelle  and  propellers  occurred  in  addition  to  the  antenna 
loss. 

A  significant  conclusion  froa  the  30*  bank  turn  exaaple  is  the  laportance 
of  tracking  all  visible  satellites  rather  than  a  subset  of  four  SVs.  Despite 
the  momentary  drop  of  SV5,  there  were  always  at  least  four  satellites  in  track 
during  the  turn.  As  a  result,  the  position  solution  was  aaintained 
continuously  throughout  the  turn. 

The  results  froa  the  flight  tests  using  the  preliainary  position  software 
can  be  summarized.  First,  the  accuracy  of  the  GPS  navigator  exceeded  the 
requirements  of  AC90-45A  for  two-dimensional  area  navigation  in  the  enroute, 
terminal  and  non-precision  approach  modes.  Second,  the  system  aaintained  the 
position  solution  during  30*  bank  turns.  Finally,  it  is  seen  that  the  system 
can  maintain  continuous  position  updates  during  momentary  outages  due  to  a 
satellite  fade  if  adequate  satellite  coverage  is  provided. 

4.3. 1.2  Engineering  Tests  with  Final  Position  Software 

The  position  software  was  frozen  at  version  2.1  at  the  beginning  of 
December,  1982.  The  remaining  engineering  flight  tests  were  conducted  with 
this  software  version.  The  ala  of  these  tests  were  to:  1)  measure  the 
accuracy  improvement  in  the  version  2.1  software,  2)  verify  performance  during 
takeoffs  and  landings,  3)  examine  the  behavior  of  the  a-0  tracker  and  4) 
determine  the  possible  effect  of  multipath  on  system  performance. 


4. 3. 1.2.1  Accuracy  Tests 

To  illustrate  the  accuracy  performance  achieved  using  the  final  version 
of  the  position  software,  the  Hanscom  flight  shown  in  Fig.  4-35  was  selected. 
The  flight  was  one  of  several  touch-and-go  landings  performed  on 
12  December  1982  at  runway  29  and  the  results  are  summarized  in  Table  4-11. 

For  this  flight  segment  lasting  about  ten  ainutes,  the  ras  horizontal  error 
was  104  feet  and  the  95Z  horizontal  error  was  175  feet.  It  should  also  be 
noted  that  the  95Z  vertical  error  in  this  case  was  215  feet,  which  aarginally 
meets  the  vertical  navigation  requirements  of  AC90-45A  shown  in  Table  2-3. 

The  horizontal  and  vertical  errors  for  the  12  December  1982  flight  segment  are 
also  shown  in  Figs.  4-36a  and  4-36b  for  the  position  fixes  and  the  tracker 
estlaates,  respectively.  The  two  large  horizontal  errors  at  about  150  sec 
into  the  run  are  due  to  aisses  in  the  HDDSEF  surveillance  rather  than  GPS 
errors;  these  artifacts  do  not  significantly  affect  the  error  statistics. 


POSITION  ACCURACY 

TABLE  4-11 

STATISTICS,  21  DECEMBER 

1982  TEST 

Position  Fixes: 

Standard 

Error 

Mean  Deviation 

RMS 

Horizontal 

34 

57 

66 

112  ft 

Vertical 

65 

49 

81 

162  ft 

Tracker  Estimates: 

Error 

Mean 

Standard 

Deviation 

RMS 

95% 

175  ft 
215  ft 


TEST 


Figure  4-37  shows  the  pseudo-range  residuals  for  the  21  December  1982 
flight*  The  pseudo-range  residuals  denote  the  difference  between  the  measured 
pseudo-range  and  the  calculated  range  based  on  the  a -8  tracker  position 
estimate.  Note  that  the  magnitudes  of  the  residuals  Increase  for  roughly  100 
seconds  starting  at  about  400  seconds  into  the  run.  This  time  interval 
corresponds  to  the  sharp  turn  during  the  flight  near  the  Shaker  Hills  waypoint 
(SKR  NDB). 

The  pseudo-range  residuals  increase  during  turns  because  of  the  behavior 
of  the  a-0  tracker.  The  tracker  projects  the  last  position  estimate  forward 
to  the  next  position  fix  time,  assuming  that  the  aircraft  is  moving  in  a 
straight  line  at  constant  velocity.  Because  acceleration  is  not  modelled,  the 
tracker  position  estimate  error  increases  during  a  turn.  As  a  result,  the 
pseudo-range  residuals  also  increase,  forcing  the  position  solution  algorithm 
to  apply  larger  corrections  to  the  position  estimate  in  order  to  produce  the 
position  fix. 

4.3. 1.2.2  Tracker  Performance 

As  a  further  illustration  of  the  effect  of  the  a-£  tracker  on  the 
pseudo-range  residuals  during  turns  consider  the  flight  segment  of  Fig.  4-38. 
In  this  test,  a  series  of  five  360s  turns  was  made  at  15s  to  20s  bank  angle. 
Figure  4-39  shows  the  pseudo-range  residuals  for  this  test.  It  is  readily 
apparent  from  the  plot  that  the  five-turn  series  causes  a  variation  in  the 
pseudo-range  residuals  of  about  +  100  feet.  It  may  be  further  noted  that  the 
residuals  for  SV8  vary  in  the  opposite  sense  from  the  other  satellites;  this 
behavior  stems  from  the  fact  that  SV8  is  west  of  the  aircraft  while  the  rest 
of  the  satellites  lie  to  the  east. 

4.3. 1.2.3  Multipath  Tests 

Tests  were  also  conducted  to  determine  if  multipath  effects  could  be 
observed.  In  one  multipath  test  the  aircraft  was  flown  over  the  ocean  on  a 
generally  south-westerly  course  while  descending  from  1500'  to  500'  in 
altitude  and  then  climbing  back  to  2500'.  The  satellite  elevation  angles 
during  this  test  ranged  from  25s  to  55*.  The  flight  profile  was  selected  in 
order  to  maximize  multipath  effects,  which  should  vary  according  to  altitude. 
The  sea  state  was  calm,  also  promoting  maximum  multipath  magnitude.  The 
ground  path  and  altitude  profile  for  the  test  are  shown  in  Figs.  4-40  and 
4-41,  respectively.  The  satellite  visibilities  are  shown  in  Fig.  4-42a. 

The  pseudo-range  residuals  for  the  test  are  shown  in  Fig.  4-42b.  It  can 
be  seen  that  the  residuals  do  not  change  in  character  throughout  the  test. 
There  is  a  bias  in  the  SV9  residuals  of  about  50  feet,  but  this  bias  does  not 
change  with  altitude. 


It  was  concluded  from  this  test  that  multipath  effects  were  not 
significant  for  satellites  of  25s  elevation  angle  or  greater.  There  were  no 
effects  discovered  in  any  of  the  flight  tests  that  could  be  attributed  to 
multipath. 
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FIG.  4-39.  PSEUDO  RANGE  RESIDUALS  FOR  FIVE  360*  TURNS 


Based  on  the  engineering  tests,  several  conclusions  were  reached.  First, 
the  accuracy  of  the  GPS  systen  exceeds  the  requirements  of  AC90-45A  for 
two-dimensional  area  navigation  in  the  enroute,  terminal  and  non-precision 
approach  modes.  As  expected,  the  system  accuracy  was  not  sufficient  to  meet 
all  of  the  requirements  of  AC90-45A  for  precision  (three-dimensional)  ILS-type 
landings,  although  the  performance  was  surprisingly  close  to  meeting  some  of 
those  requirements.  The  GPS  accuracy  also  met  the  projected  requirements  for 
non-precision  approach  as  stated  in  the  Federal  Radio  Navigation  Plan. 

Finally,  the  system  performance  exceeded  the  project  goals  for  level  and  30s 
bank  turns. 

The  second  conclusion  from  the  engineering  tests  was  that  the  system  is 
not  unacceptably  affected  by  aircraft  dynamics.  The  ability  of  the  system  to 
maintain  position  updates  throughout  30°  bank  angle  turns  was  demonstrated. 

It  was  also  shown  that  the  system  maintained  continuous  position  updates 
during  momentary  satellite  outages  during  a  turn. 

A  third  conclusion  was  that  the  a-B  tracker  performance  is  adequate  even 
though  it  does  not  attempt  to  model  vehicle  acceleration.  Given  the  low 
airspeed  and  limited  dynamics  of  general  aviation  aircraft,  the  use  of  the 
linear  tracker  appears  justified.  The  test  results  show  that  the  system 
accuracy  remains  well  within  the  horizontal  accuracy  requirements  of 
AC90-45A. 

A  final  conclusion  from  the  engineering  tests  was  that  multipath 
interference  did  not  appear  to  have  a  significant  effect  on  system  performance 
Despite  specific  tests  designed  to  produce  multipath  effects,  no  pseudo-range 
errors  occurred  which  could  be  attributed  to  multipath. 

4.3.2  Operational  Flight  Tests 

Upon  completion  of  the  engineering  flight  tests,  a  series  of  operational 
flights  were  conducted  to  test  the  performance  of  the  system  as  an  area 
navigator.  These  operational  tests  were  conducted  in  three  different  areas 
representing  various  airport  types  likely  to  be  encountered  by  general 
aviation  aircraft.  The  areas  were: 

•  Burlington  International  Airport  -  Burlington,  Vt. 

•  Logan  International  Airport  -  Boston,  MA. 

•  Hans com  Field  -  Lexington  -  Manchester  Airport  -  Manchester,  NH. 

The  tests  at  the  Burlington,  Vt.  airport  were  designed  to  determine  if  a  GPS 
navigator  can  operate  in  an  area  surrounded  by  mountainous  terrain.  The  tests 
at  Logan  International  were  made  to  verify  that  GPS  can  be  operated  at  an 
urban  airport  adjacent  to  large  buildings.  The  flights  in  the  Hanscom 
Field/Manchester  Airport  area  were  designed  to  verify  the  correct  performance 
of  the  GPS  navigation  system  in  typical  general  aviation  operations. 


4.3.2. 1  Burlington,  Vt.  Tests 

The  approach  patterns  to  Burlington  International  airport  are  shown  in 
Fig.  4-43.  There  are  two  non-precision  approaches  to  the  airport:  an  NDB 
approach  to  Runway  IS  and  a  VOR  approach  to  Runway  1.  The  airport  is  located 
between  two  mountain  ranges,  the  Green  Mountains  at  about  20  ailes  to  the  east 
and  the  Adirondack  Mountains  at  about  25  ailes  to  the  southwest. 

For  the  Burlington  operational  tests,  the  GPS  test  aircraft  was  flown 
froa  Hans cob  Field  in  Massachusetts  to  Burlington  International  on 
27  January  1983.  A  series  of  non-precision  approaches  were  aade  using  both 
the  VOR  and  NDB  approach  procedures.  Both  the  BTV  VOR  and  the  BT  NDB  were 
entered  as  waypoints  in  the  GPS  navigation  software  for  this  purpose. 

4. 3. 2. 1.1  Enroute  Performance 


During  the  first  portion  of  the  test  flight,  the  aircraft  was  enroute  to 
the  Burlington  VOR  froa  Hanscoa  Field  as  shown  in  Fig.  4-44.  At  the  beginning 
of  run  2a  the  aircraft  was  70  nautical  ailes  away  froa  the  VOR  on  the  151° 
radial.  At  the  altitude  the  aircraft  was  flying  (4250  feet),  the  VOR-DME  was 
unusable  for  ranges  greater  than  30  na  on  this  radial. 

Figure  4-45  shows  the  altitude  and  ground  speed  for  run  2a  as  estimated 
by  the  position  software.  Fig.  4-46  shows  the  GPS  satellite  azimuth  and 
elevation  during  the  run.  Note  that  only  four  satellites  were  visible  during 
this  period.  The  satellite  carrier-to-noise-density  values  are  shown  in 
Fig.  4-47.  Note  at  the  beginning  of  the  run  that  SV9  has  a  C/N0  about  5  dB 
lower  than  the  other  SVs.  At  this  tlae,  SV9  was  at  8*  elevation  angle,  but 
still  had  a  substantial  loaa-of-lock  aargin  of  about  12  dB. 

The  pseudo-range  residuals  for  2a  are  shown  in  Fig.  4-48.  The  residuals 
are  seen  to  reaaln  small,  typically  less  than  +  50  feet,  over  most  of  the  run. 
There  are,  however,  three  Instances  where  the  residuals  became  large:  once 
for  SV9  and  twice  for  SV4.  Analysis  of  the  data  showed  that  these 
perturbations  occured  when  the  position  software  was  unable  to  collect  a  valid 
pseudo-range  aeasureaent  froa  the  receiver.  Due  to  an  apparent  software  bug, 
the  failure  to  collect  a  valid  aeasureaent  caused  the  next  measurement  to  be 
laproperly  tlae-tagged.  The  improper  tlae-tagging,  in  turn,  caused  the  range 
to  the  satellite  to  be  calculated  incorrectly.  The  Incorrect  range 
calculation  then  caused  the  pseudo-range  residual  to  be  in  a  error  by  about 
500  feet. 

It  is  significant  that  a  pseudo-range  error  of  this  magnitude  has  only  a 
saall  effect  on  the  navigation  accuracy;  this  is  due  to  two  reasons.  First, 
the  error  is  short-tern,  lasting  for  a  maximum  of  two  navigation  cycles. 
Second,  the  batch  least-squares  aethod  alnlaises  the  effect  of  large  residuals 
on  the  position  fix.  In  this  case,  the  500  foot  residual  error  caused  a 
position  fix  error  of  only  about  100  feet.  A  position  estlaate  error  of  this 
magnitude  is  not  dlscernable  on  the  pilot  course  deviation  display.  It  should 
also  be  noted  that  pseudo-range  reasonableness  check  in  the  position  software 
rejects  any  residual  greater  than  1000  feet,  preventing  a  bad  position  fix 
froa  being  generated. 
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The  results  of  this  test  show  the  ability  of  the  GPS  system  to  perform 
enroute  navigation  in  mountainous  terrain  with  a  minimum  set  of  four 
satellites.  It  is  of  particular  interest  that  the  GPS  was  able  to 
significantly  outperform  VOR  in  this  mountainous  area.  GPS  was  able  to 
provide  equivalent  navigation  service  70  nmi  away  from  the  VOR  while  the  VOR 
Itself  is  unusable  at  ranges  greater  than  30  nmi.  This  performance  advantage 
of  GPS  results  from  the  fact  that  the  GPS  satellites  are  above  the  horizon, 
while  VOR  is  ground-based  and  much  more  susceptible  to  terrain  blockage. 

4.3.2. 1.2  VOR  Approach 

Runs  2b  and  3b  consist  of  non-precision  VOR  approaches  to  the  Burlington 
airport.  Figure  4-49  shows  the  ground  track  for  run  2b.  The  aircraft  first 
intercepted  the  VOR  on  the  151s  radial,  then  proceeded  outbound  on  the  216° 
radial  and  performed  a  procedure  turn.  After  intercepting  the  VOR  on  the  36° 
radial,  the  aircraft  was  then  flown  to  the  field  for  a  low  approach  at  runway 
33.  Some  gaps  will  be  noted  in  the  ground  track;  the  gaps  are  due  to  data 
recording  problems  rather  than  Interruptions  in  the  position  solution. 

Figure  4-50  shows  the  altitude  and  ground  speed  profiles  for  run  2b. 
Figure  4-51  shows  the  satellite  visibilities  for  the  run.  Although  five 
satellites  are  shown,  this  run  was  still  navigating  with  four  satellites  (SVs 
4,6,8  and  9).  The  satellite  C/N0  values  are  shown  in  Fig.  4-52.  It  is  seen 
that  the  satellite  signal  strength  drops  by  as  much  as  10  dB  in  some  turns, 
but  always  remains  above  the  33  dB-Hz  loss-of-lock  threshold.  Figure  4-53 
shows  the  pseudo-range  residuals;  note  the  increases  in  the  residuals 
corresponding  to  the  turns.  There  are  also  some  large  residuals  at  the  end  of 
the  run  for  SV4;  these  jumps  result  from  the  software  bug  previously 
described. 

The  ground  track  for  run  3b  is  shown  in  Fig.  4-54.  For  this  run,  the 
aircraft  intercepted  the  VOR  then  flew  the  216°  outbound  radial  and  made  the 
procedure  turn  as  before.  After  Intercepting  the  VOR,  an  OBS  setting  of  36s 
from  the  waypoint  was  flown.  It  was  found,  however,  that  this  resulted  in  a 
42*  course  from  the  waypoint,  taking  the  aircraft  to  the  east  of  the  end  of 
runway  33.  Because  of  air  traffic  control  considerations,  a  right  downwind 
approach  was  then  flown  to  runway  15  and  a  landing  made  for  refueling. 

Post-flight  analysis  revealed  an  error  in  the  navigation  software  which 
had  the  result  of  biasing  the  OBS  setting  by  6°,  explaining  the  Incorrect 
course  of  42°  commanded  by  the  navigation  software.  This  error  was  corrected 
for  subsequent  operational  flight  tests. 


Figure  4-55  shows  the  altitude  and  ground  speed  profiles  for  run  3b. 
Note  that  during  the  last  portion  of  the  run,  the  aircraft  is  on  the  ground 
and  taxiing  around  the  airport.  The  satellite  visibilities  are  shown  in 
Fig.  4-56.  All  five  satellites  were  in  track  during  this  run,  including  SV5 
which  was  initially  at  10*  elevation  angle.  The  satellite  C/N0  values  are 
shown  in  Fig.  4-57.  Note  that  SV5,  the  low  elevation  angle  satellite,  was 
dropped  on  two  occasions  during  the  run.  These  occasions  correspond  to  the 
procedure  turn  and  the  base  leg  turn,  at  which  times  the  aircraft  was  banked 
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FIG.  4-49.  BURLINGTON.  VT  VOR  APPROACH.  27  JAN  83 
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FIG.  4-54.  BURLINGTON,  VT  VOR  APPROACH.  RUN  3b,  27  JAN  83 


away  from  SV5.  The  position  solution  Is  maintained  continuously  in  these 
cases ,  however,  since  the  other  four  satellites  remained  In  track.  This 
example  again  Illustrates  the  importance  of  tracking  all  visible  satellites 
rather  than  the  minimum  set  of  four. 

Figure  4-58  shows  the  pseudo-range  residuals  for  run  3b.  It  can  be  seen 
from  the  figure  that  there  are  a  number  of  jumps  in  the  pseudo-range 
residuals.  Except  for  the  two  cases  in  which  SV5  suffered  a  loss  of  lock,  the 
jumps  are  all  due  to  the  position  software  problem  discussed  earlier.  In  this 
run,  the  position  software  appears  to  occasionally  fail  to  collect  a  valid 
pseudo-range  measurement  before  it  is  overwritten  by  the  receiver  with  the 
next  measurement.  This  problem  is  symptomatic  of  a  lack  of  available  CPU  time 
in  the  position  processor.  The  lack  of  CPU  time  is  attributable  in  part  to 
the  large  amount  of  engineering  data  being  recorded  by  the  system.  The  data 
recording  function  introduces  an  overhead  which  would  not  be  present  in  a 
production  system. 

It  can  also  be  noted  from  Fig.  4-58  that  there  is  a  substantial  bias 
(-75  feet)  on  the  SV4  residuals  and  a  smaller  bias  in  the  opposite  direction 
on  SV6.  The  source  of  this  bias  is  not  known,  but  it  may  be  due  to  aging 
ephemeric  data  from  SV4.  Unlike  the  other  satellites,  SV4  was  operating  on  a 
crystal  oscillator  because  its  atomic  clocks  had  failed.  Consequently,  SV4's 
clock  drift  is  not  well  modelled  by  the  ephemeric  data  from  the  satellite. 

The  bias  is  clearly  not  due  to  multipath  because  it  does  not  vary  with 
altitude,  remaining  constant  throughout  the  run. 

The  results  of  runs  2b  and  3b  show  that  GPS  can  be  used  to  simulate  a  VOR 
approach  in  a  mountainous  area.  Furthermore,  despite  a  position  software  bug 
that  introduced  pseudo-range  jumps,  the  system  maintained  reliable  operation 
throughout  standard-rate  turns  even  during  the  temporary  loss  of  a  low-angle 
satellite.  Also,  the  system  kept  a  low  (10°)  elevation  angle  satellite  in 
track  during  level  flight. 

4. 3. 2. 1.3  NDB  Approach 

The  ground  track  for  two  NDB  approaches  to  runway  15  (run  4)  is  shown  in 
Fig.  4-59.  For  this  run,  the  aircraft  took  off  from  runway  15  and  proceeded 
outbound  on  the  326°  radial  from  the  NDB.  After  a  procedure  turn,  the 
aircraft  made  a  non-precision  approach  on  the  146°  radial.  A  missed  approach 
was  executed  on  advice  from  air-traffic  control  and  a  second  approach 
initiated.  At  this  point,  the  run  terminated  due  to  a  momentary  aircraft 
electrical  power  interruption. 


The  altitude  and  ground  speed  profiles  for  run  4  are  shown  in  Fig.  4-60. 
The  satellite  visibilities  are  shown  in  Fig.  4-61.  Figure  4-62  shows  the 
measured  satellite  C/N0  values  for  the  turn.  Note  the  large  number  of  losses 
of  lock  for  SV6  during  the  run.  Examining  the  satellite  elevation,  aircraft 
altitude  and  topography  of  the  region,  it  is  clear  that  these  dropouts  are  not 
due  to  terrain  blockage.  Examination  of  the  data  shows  that  SV6  was  dropped 
once  temporarily  during  the  takeoff  due  an  aircraft  power  surge.  Lock  was 
reestablished  for  the  satellite  in  the  first  SV6  dwell  slot  in  the  navigation 
cycle  but  not  in  the  second.  As  a  result,  the  satellite  C/N0  estimates  for 
the  second  dwell  slot  are  not  reliable,  but  were  lnadvertenly  plotted. 
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FIG.  4-69.  BURLINGTON,  VT  NDB  APPROACH,  RUN  4,  27  JAN  83 
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ALTITUDE  AND  GROUND  SPEED,  RUN  4,  27  JAN  83 


SATELLITE  C/NO.  VALUES,  RUN  4.  27  JAN  83 


It  should  be  noted  that  the  verlson  2.1  position  software  does  not 
contain  a  provision  for  replacing  failed  satellites  in  the  navigation  mode 
dwell  cycle.  This  provision  was  built  into  version  3.0  software  which  was  not 
completed  in  time  to  support  flight  testing.  Had  the  satellite  replacement 
feature  been  Installed,  the  position  software  would  have  been  able  to 
reestablish  lock  in  the  second  SV6  dwell  slot.  It  should  be  noted  that  the 
satellite  replacement  feature,  while  not  flight  tested,  was  validated  for 
version  3.0  in  laboratory  tests. 

The  pseudo-range  residuals  for  run  4  are  shown  in  Fig.  4-63.  The 
residuals  are  similar  to  those  for  run  3.  Again  note  that  SV4  had  a  negative 
bias  and  SV6  had  a  positive  bias.  The  SV6  positive  bias  is  probably  due  to 
the  batch  least-squares  position  fixing  algorithm  attempting  to  minimize  the 
negative  bias  on  SV4  due  to  its  less  accurate  clock. 

The  results  for  run  4  are  seen  to  be  similar  for  the  VOR  approaches.  It 
is  seen  that  the  GPS  system  functioned  well  on  the  NDB  approach  in  this 
mountainous  region. 


4. 3. 2. 2  Logan  International  Tests 

Operational  GPS  test  flights  were  made  into  Logan  International  airport 
at  Boston,  HA  on  25  January  1983  and  28  January  1983.  The  number  of  flights 
into  Logan  was  limited  by  difficulty  in  obtaining  ATC  clearances  into  the 
airport  during  the  satellite  visibility  period,  which  coincided  with  time  of 
high  air  traffic  density.  However,  enough  experience  was  gained  at  Logan  to 
verify  successful  operation  of  the  GPS  system  at  this  large,  busy  urban 
airport. 

Figure  4-64  shows  the  ground  track  recorded  by  the  GPS  system  during  run 
2  on  25  January  1983.  The  test  aircraft  approached  the  LYNDY  waypoint  (a 
non-directlonal  beacon  site)  from  the  northeast.  After  intercepting  the 
waypoint,  the  aircraft  began  a  non-precison  approach  to  runway  22.  Upon 
landing,  the  aircraft  taxlied  to  the  end  of  runway  22R  and  took  off.  After 
intercepting  the  MILTT  waypoint  (also  an  NDB  site),  the  aircraft  turned  to  the 
northwest  and  departed  the  Logan  area.  Throughout  the  run,  the  GPS  system 
maintained  all  five  satellites  in  track  and  supplied  continuous  position 
updates.  (Note:  the  gaps  in  the  ground  track  shown  in  Fig.  4-64  are  due  to 
data  recording  dropouts). 

Figure  4-65  shows  the  altitude  and  ground  speed  recorded  by  the  GPS 
system  during  the  Logan  run.  As  seen  in  the  figure,  the  middle  part  of  the 
run  was  spent  taxiing  around  the  airport.  During  the  time  the  aircraft  was 
taxiing,  there  were  numerous  airport  buildings  and  large  aircraft  in  the 
immediate  vicinity.  Figure  4-66  shows  the  satellite  visibilities  for  the  run. 
Note  that  four  of  the  satellites  are  to  the  west  of  the  airport,  towards  the 
downtown  Boston  area,  and  the  remaining  satellite  is  to  the  east,  towards  the 
ocean.  Figure  4-67  shows  the  dilution-of -precision  values  for  the  run  an 
calculated  in  post-processing.  As  seen  in  the  figure,  the  GDOP  varied  from 
about  6.3  to  6.7  during  the  run. 
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The  satellite  C/N0  values  are  shown  In  Fig.  4-68.  All  of  the  satellite 
C/N0  values  remained  well  above  the  33  dB-Hz  loss-of-lock  threshold,  although 
SV6  dropped  substantially  (about  10  dB)  during  the  turn  after  takeoff  from 
runway  22  R.  The  drop  in  0/No  is  not  surprising  since  the  aircraft  banked 
away  from  SV6  during  the  turn. 

It  is  significant  that  the  GPS  system  kept  the  satellites  in  track 
continuously  while  taxiing  despite  the  buildings  and  large  aircraft  near  to 
the  GPS  aircraft.  Although  blockage  of  the  satellite  signal  is  probable  when 
the  aircraft  is  immediately  adjacent  to  a  building,  significant  blockage  did 
not  occur  on  the  taxiways  or  runways. 

The  pseudo-range  residuals  for  the  25  Janaury  1983  run  are  shown  in 
Fig.  4-69.  The  residuals  showed  no  jumps  of  the  kind  seen  earlier,  indicating 
that  the  position  software  was  always  able  to  collect  pseudo-range 
measurements  from  the  receiver.  Note  the  positive  bias  (+75  feet  on  SV9  and 
negative  bias  (-50  feet)  on  SV5.  As  explained  earlier,  these  biases  appear  to 
be  due  to  old  ephermeric  data  transmitted  by  the  satellites. 

A  second  test  flight  was  made  at  Logan  International  on  28  January  1983. 
Figure  4-70  shows  the  ground  track  for  a  low-approach  and  departure  from 
runway  4L.  The  altitude  and  ground  speed  profiles  for  the  flight  are  shown  in 
Fig.  4-71.  After  departing  the  runway,  the  aircraft  made  a  turn  to  the 
northwest  prior  to  intercepting  the  LTNDY  waypoint  on  instructions  from  air 
traffic  control.  The  aircraft  then  departed  the  Logan  area  for  Hans com 
Field. 

The  satellite  C/N0  values  for  the  28  January  1983  Logan  run  are  shown  in 
Fig.  4-72.  The  run  was  made  late  in  the  satellite  visibility  window,  so  that 
only  four  satellites  were  available  for  navigation.  Nonetheless,  the  GPS 
system  was  able  to  provide  continuous  navigation  updates  throughout  the  run. 
The  pseudo-range  residuals  for  the  run  are  shown  in  Fig.  4-73.  Note  the 
positive  change  la  the  SV8  residual  and  negative  changes  in  the  SV5  and  SV9 
residuals  corresponding  to  the  initial  turn  in  the  run.  Again,  this  behavior 
is  explained  by  the  satellite  positions:  SV8  was  to  the  east  while  SV5  and 
S "ere  to  the  west.  SV4  was  more  to  the  north  and  the  aircraft  velocity  in 
that  direction  changed  little. 

The  Logan  runs  demonstrated  that  GPS  can  be  used  for  non-precision 
approach  into  a  large  urban  airport.  It  was  also  found  that  the  GPS  system 
was  able  to  maintain  navigation  service  while  on  the  taxiways  and  runways 
under  conditions  that  generally  render  VOR/DME  systems  inoperable.  It  should 
be  cautioned  that  these  observations  apply  to  the  GPS  system  while  in 
navigation  mode;  the  lock  threshold  is  2  dB  higher  (less  sensitive)  for  the 
acquisition  mode.  Also,  some  satellite  signals  are  likely  to  be  blocked  when 
the  aircraft  is  parked  next  to  a  large  building  such  as  a  passenger  terminal. 
Nonetheless,  the  satellite  C/N0  a*r8la*  appear  quite  acceptable  once  the 
aircraft  is  out  on  the  taxiways. 
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4. 3. 2. 3  Hanscon  Field /Manchester  Airport  Tests 


4. 3. 2. 3.1  25  January  1983  Teat 

The  final  set  of  operational  tests  were  nade  in  the  vicinity  of  Hanscon 
Field  at  Lexington,  MA  and  Manchester  Airport  at  Manchester,  NH.  The  purpose 
of  these  tests  was  to  verify  the  use  of  GPS  as  a  navigation  system  in  typical 
general  aviation  operations.  Figure  4-74  shows  the  GPS-generated  ground  track 
of  such  a  flight  into  Manchester  Airport  on  25  January  1983.  In  this  case  the 
aircraft  was  flown  north  until  the  337*  radial  to  the  MHT  VOR  waypoint  was 
intercepted.  The  aircraft  was  then  flown  on  the  radial,  intercepting  the 
waypoint  and  descending  for  a  landing  at  runway  35. 

The  altitude  and  ground  speed  profiles  for  the  flight  are  shown  in 
Fig.  4-75.  As  seen  in  the  figure,  a  substantial  disturbance  occured  in  the 
navigation  updates  just  prior  to  landing.  Analysis  of  the  test  data  showed 
that  this  disturbance  was  caused  by  an  incorrect  position  fix  calculation. 

The  calculation  was  Incorrect  because  the  position  software  failed  to  collect 
valid  measurements  from  the  receiver  in  both  SV4  dwells  and  a  pseudo-range 
from  the  previous  navigation  cycle  was  included  in  the  position  fix 
calculation  due  to  a  position  software  error.  The  system  corrected  the 
position  error  on  the  next  cycle,  2.2  seconds  later. 

It  should  also  be  noted  that  the  magnitude  of  the  change  in  the  position 
fix  was  greater  than  a  thousand  feet,  so  that  the  fix  should  have  been 
declared  Invalid  by  the  position  software.  However,  the  limits  on  position 
fix  changes  were  unrealistically  high  in  the  Version  2.1  software  used  for  the 
test,  so  that  the  fix  was  passed  on  to  the  a-p  tracker.  Despite  this  fact, 
the  disturbance  was  sufficiently  damped  out  such  that  it  was  not  operationally 
significant. 

The  satellite  visibilities  for  the  25  January  1983  run  are  shown  in 
Fig.  4-76.  Note  that  SV6  was  successfully  tracked  down  to  4°  elevation  angle, 
below  the  nominal  5*  elevation  angle  mask.  Even  more  surprising,  the  aircraft 
was  on  the  ground  at  the  time.  Figure  4-77  shows  the  satellite  C/N0  values 
for  run.  Note  that  SV6  is  dropped  for  a  short  time,  just  before  landing,  but 
1 8  brought  back  into  track.  Shortly  afterward,  SVs  4,5,8  and  9  are  also 
dropped  but  tracking  is  immediately  recovered. 

Figure  4-78  shows  the  pseudo-range  residuals.  The  jump  in  the  SV4 
residual  was  caused  by  the  incorrect  position  fix,  as  previously  discussed. 

The  other  two  jumps  occur  in  SV8  and  are  caused  by  the  temporary  loss  of 
tracking  for  that  satellite. 


4. 3. 2. 3. 2  4  February  1983  Test 

For  subsequent  operational  tests  in  the  Hanscom  Field /Manchester  Airport 
area,  a  standard  test  scenario  was  developed  as  shown  in  Fig.  4-79.  For  this 
scenario,  the  aircraft  departs  from  Hanscon  Field,  runway  11,  and  is  flown  to 
the  Shaker  Hills  NDB  waypoint  (SKR  4).  After  intercepting  the  waypoint,  the 
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aircraft  then  Intercepts  the  356°  radial  to  the  Manchester  VOR  waypoint 
(MHT  6).  After  overflying  the  Manchester  waypoint,  the  aircraft  Is  turned  to 
157°  outbound,  performs  a  procedure  turn  and  turns  inbound  on  the  337*  radial. 
Upon  intercepting  the  waypoint  again,  the  aircraft  begins  the  descent  to 
runway  35.  A  missed  approach  is  then  performed  with  the  Manchester  waypoint 
as  the  Missed  Approach  Point  (MAP). 

After  overflying  the  Manchester  waypoint  again,  the  aircraft  departs  on 
the  221*  radial  for  the  LOBBY  intersection  waypoint  (LOBBY  2).  The  aircraft 
enters  a  holding  pattern  at  LOBBY,  then  departs  on  the  113*  radial  for  the 
Bedford  NDB  waypoint  (BE  1).  The  descent  to  runway  11  at  Hanscom  Field  begins 
after  the  Bedford  waypoint  is  overflown.  The  aircraft  then  lands  at  Hanscom 
or  performs  a  low-approach  to  start  another  circuit  of  the  course. 

The  ground  track  for  Run  la  on  4  February  1983  is  shown  in  Fig.  4-80. 

The  test  aircraft  was  under  continuous  MODSEF  surveillance  commencing  just 
before  the  aircraft  reached  the  SKR  waypoint  and  ending  as  the  aircraft  left 
the  LOBBY  waypoint.  As  can  be  seen  from  the  figure,  the  aircraft  was 
successfully  navigated  to  the  selected  waypoints  using  the  GPS  system.  On  the 
Shaker  Hills  to  Manchester  leg,  the  pilot  maintained  the  course  indicated  by 
the  system  fairly  closely,  but  departed  from  the  GPS  course  on  the  Manchester 
to  LOBBY  leg  due  to  ATC  traffic  advisories.  The  altitude  and  ground  speed 
profiles  for  the  run  are  shown  in  Fig.  4-81. 

The  satellite  visibilities  for  Run  la  are  shown  in  Fig.  4-82  and  the 
dilution-of -precision  calculations  are  given  in  Fig.  4-83.  Figure  4-84  shows 
the  satellite  C/N0  values  during  the  run  and  Fig.  4-85  shows  the  pseudo-range 
residuals. 

Since  the  test  aircraft  was  under  MODSEF  surveillance  during  the  run  the 
system  accuracy  could  be  determined  by  post-flight  processing.  Figure  4-86a 
shows  the  horizontal  and  vertical  position  fixing  error  as  a  function  of  time. 
There  are  several  large  horizontal  errors  shown  in  the  figure;  analysis  of  the 
data  shows  that  the  large  values  are  instances  in  which  the  MODSEF  facility 
failed  to  track  the  aircraft.  These  spurious  errors  do  not  substantially 
affect  the  error  statistics.  The  horizontal  and  vertical  error  in  the  GPS 
position  estimate  (tracker  output)  is  shown  in  Fig.  4-86b. 

Table  4-12  shows  the  position  accuracy  statistics  for  Run  la.  The 
horizontal  accuracy  of  the  tracker  estimates  is  seen  to  be  well  within  the 
system  performance  requirements  and  goals. 

In  order  to  examine  the  behavior  of  the  system  during  a  turn,  the  holding 
pattern  section  of  Run  la  was  expanded  as  shown  in  Figs.  4-87  and  4-88.  The 
MODSEF  radar  track  is  plotted  with  the  GPS  position  fixes  in  Fig.  4-87  and 
with  the  tracker  estimates  in  Fig.  4-88.  The  horizontal  and  vertical  errors 
for  the  position  fixes  and  tracker  estimates  are  shown  in  Figs.  4-89  and  4-90, 
respectively.  As  seen  in  Fig.  4-87  and  4-89  the  horizontal  position  fix  error 
does  not  build  up  during  the  turns,  and  the  only  major  discrepancies  between 
the  GPS  position  fixes  and  the  radar  track  appear  to  be  MODSEF  tracking 
errors.  Note  in  particular  that  during  Turn  3  the  horizontal  position  fix 
error  remains  less  than  150  feet  throughout  the  turn. 
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FIG.  4-87.  HOLDING  PATTERN,  4  FEB  83.  RUN  1«:  POSITION  FIXES 
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Figures  4-88  and  4-90  show  the  performance  of  the  tracker  in  the  holding 
pattern.  It  is  seen  that  a  considerable  bias  error  occurs  in  the  tracker 
estimates  during  some  of  the  turn  segments.  This  bias  error  is  especially 
noticeable  in  Turn  3  where  the  horizontal  error  Increases  to  as  much  as 
300  feet.  A  possible  explanation  for  the  increased  bias  error  in  Turn  3  can 
be  obtained  from  Fig.  4-91  which  shows  the  altitude  and  ground  speed  profile 
for  the  holding  pattern.  It  is  seen  that  the  ground  speed  was  about  300  feet 
per  second  (180  kts)  at  the  beginning  of  Turn  3  but  declined  to  200  fps  by  the 
end  of  the  turn  60  seconds  later.  Because  this  deceleration  is  unmodelled, 
the  tracker  tended  to  overestimate  the  distance  traveled  between  fixes  at  the 
beginning  of  the  turn.  As  the  ground  speed  declined,  the  errors  in  Turn  3 
also  declined.  By  contrast,  the  ground  speed  was  low  at  the  start  of  Turn  4 
but  increased  during  the  turn;  consequently,  the  horizontal  errors  tended  to 
Increase  during  the  turn. 

The  position  error  statistics  for  the  holding  pattern  are  shown  in 
Table  4-13.  The  horizontal  error  of  the  tracker  estimates  is  clearly  well 
within  the  system  requirements  and  goals. 

Run  lb  on  4  February  1982  was  similar  to  Run  la,  except  that  a  low 
approach  was  made  at  Manchester  airport  and  the  aircraft  was  flown  back  to  the 
Bedford  NDB  (BE  1)  waypoint  instead  of  the  LOBBY  waypoint.  The  accuracy 
statistics  were  similar  to  Run  la. 

4. 3. 2. 3. 3  9  February  1983  Test 

A  second  operational  test  of  the  type  conducted  on  4  February  1983  was 
made  on  9  February  1983.  Figure  4-92  shows  the  ground  track  for  Run  la.  For 
this  test,  the  aircraft  took  off  from  Hanscom  Field  runway  29  (instead  of 
runway  11)  due  to  prevailing  wind  conditions.  After  departing  the  runway,  the 
aircraft  was  turned  180°  toward  the  Shaker  Hills  NDB  (SKR).  The  remainder  of 
the  run  was  made  according  to  the  test  scenario  of  Fig.  4-79.  The  gap  in  the 
ground  track  after  the  turn  at  the  SKR  NDB  is  due  to  an  inadvertent 
Interruption  in  the  data  recording;  the  GPS  navigation  updates  were  continuous 
throughout  the  run.  The  altitude  and  ground  speed  profile  of  the  run  is  shown 
in  Fig.  4-93. 

The  satellite  visibilities  are  shown  in  Fig.  4-94.  Although  the 
visibilities  for  all  five  satellites  are  shown,  only  four  satellites  (SVs 
4,5,6  and  8)  were  used  for  navigation  during  the  run.  Despite  this  the  system 
maintained  continuous  navigation  updates  for  a  period  of  nearly  two  hours. 

The  dllutlon-of-precision  calculations  for  Run  la  are  shown  in  Fig.  4-95. 

Note  that  the  GDOP  for  the  run  is  higher  than  that  for  the  test  of  4  February 
1983  due  to  the  use  of  only  four  satellites. 


The  horizontal  and  vertical  position  fix  error  versus  time  is  shown  in 
Fig.  4-96.  The  correspondihg  data  for  the  tracker  estimates  is  shown  in 
Fig.  4-97.  It  should  be  noted  that  the  M0DSEF  radar  tracking  data  did  not 
start  until  time  ■  43085  sec  due  to  difficulties  with  the  radar  system;  the 
start  of  the  radar  tracking  is  indicated  in  Fig.  4-92.  The  accuracy 
statistics  for  the  run  are  summarized  in  Table  4-14.  These  statistics  compare 
favorably  to  the  4  February  data. 


TABLE  4-13 


POSITION  ACCURACY  STATISTICS  FOR  HOLDING  PATTERN 
4  FEBRUARY  1983,  RUN  1A 


Position  Fixes: 

Standard 

Error 

Mean 

Deviation 

RMS 

Horizontal 

55 

75 

93 

Vertical 

151 

73 

198 

Tracker  Estimates: 

Standard 

Error 

Mean 

Deviation 

RMS 

Horizontal 

49 

49 

70 

Vertical 

157 

98 

185 

95Z 
155  ft 
277  ft 


95Z 
253  ft 
331  ft 
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400- 


-400- 


TABLE  4-14 


POSITION  ACCURACY  STATISTICS  FOR  9  FEBRUARY  1983 
RUN  1A  (43085  -  45770) 


Position  Fixes: 


Standard 

Error 

Mean 

Deviation 

RMS 

Horizontal 

73 

112 

134 

Vertical 

68 

48 

83 

Tracker  Estimates: 


Error 


Standard 

Mean  Deviation  RMS 


At  the  end  of  run  la,  the  aircraft  had  completed  one  complete  circuit 
around  the  test  course  so  a  second  circuit  as  started  (run  lb)  as  shown  in 
Fig.  4-98.  The  aircraft  made  a  low-approach  at  500  feet  over  Hanscom  Field 
and  proceeded  to  the  Shaker  Hills  waypoint.  The  run  terminated  near  the  LOBBY 
waypoint  when  SV6  dropped  below  the  5s  elevation  angle  mask.  The  altitude  and 
ground  speed  profiles  for  run  lb  are  shown  in  Fig.  4-99. 

The  satellite  visibilities  for  run  lb  are  shown  in  Fig.  4-100.  Again 
note  that  SV9 ,  although  Included  in  the  figure,  was  not  used  for  navigation. 
Mote  also  at  the  end  of  the  run  that  SV6  had  dropped  below  the  nominal  5° 
elevation  angle  mask.  The  dilution-of-precislon  values  for  the  run  are  shown 
in  Fig.  4-101. 

Figure  4-102  shows  the  satellite  C/N0  values  for  the  latter  portion  of 
run  lb.  Note  that  the  SV6  C/NQ  value  at  the  end  of  the  run  is  above  the 
33  dB-Hz  loss-of-lock  threshold.  The  satellite  C/NQ  does  drop  below  the 
threshold  several  times  during  the  initial  turn  of  the  holding  pattern  but  the 
satellite  is  quickly  re-established  in  track  on  each  occasion. 

The  horizontal  and  vertical  position  fix  error  vs.  time  is  shown  in 
Fig.  4-103.  The  corresponding  tracker  estimate  error  is  shown  in  Fig.  4-104. 
The  position  accuracy  statistics  are  summarized  in  Table  4-15.  In  general, 
the  position  accuracy  is  seen  to  be  less  for  run  lb  than  Run  la;  however,  the 
GDOP  has  also  become  larger  in  the  latter  run  as  seen  from  Fig.  4-101. 
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ALTITUDE  AND  GROUND  SPEED,  9  FEB  83,  RUN  1b  (2  OF  2) 


223 


PRECISION  VALUES.  9  FEB  83.  RUN  1b 


600 


-300 


TABLE  4-15 


POSITION  ACCURACY  STATISTICS  FOR  9  FEBRUARY  1983 
RUN  IB  (46030  -  48580) 


Position  Fixes: 

Error 

Mean 

Standard 

Deviation 

RMS 

Horizontal 

185 

117 

219 

Vertical 

153 

83 

174 

I 


Tracker  Estimates: 


Standard 

Deviation 

RMS 

Horizontal 

191 

142 

238 

Vertical 

159 

97 

186 

95% 
322  ft 
293  ft 


95% 
378  ft 
329  ft 


227 


4.3. 2. 4  Su— ry  of  Operational  Plight  Test  Results 

The  operational  flight  teats  demonstrated  the  use  of  the  GPS  C/A  code 
navigator  In  mountainous  terrain,  at  a  large  urban  airport  and  in  typical 
general  aviation  operations.  The  position  accuracy  for  the  operational  flight 
tests  In  the  Hanscoa  Field/Manchester  Airport  is  suamarlaed  in  Table  4-16. 
These  statistics  combine  the  results  from  4  February  1983  (Run  la), 

9  February  1983  (Runs  la  and  lb)  and  10  February  1983  (Run  2)  into  a  single 
run  lasting  2.7  hours.  As  seen  from  the  table,  the  95Z  horizontal  accuracy  of 
the  tracker  estiaates  was  333  feet,  which  easily  meets  the  0.3  na  accuracy 
requirement  of  AC90-45A  for  two-dimensional  area  navigation.  This 
performance  also  essentially  aeets  the  propoaed  328  foot  (95Z)  accuracy 
requirements  of  the  Federal  Navigation  Plan  for  non-precision  approach. 

The  operational  flight  tests  also  demonstrated  the  compatibility  of  the 
GPS  systea  with  current  air  navigation  practicea  and  procedurea.  The  teat 
pilots  responded  favorably  to  the  operation  of  the  systea,  noting  in 
particular  that  the  GPS  navigator  gave  a  more  stable  cross-track  deviation 
indication  for  non-precision  approach  than  that  provided  by  a  VOR  receiver. 

It  was  also  noted  that  the  GPS  could  obtain  the  course  and  distance  to  a 
selected  waypoint  while  on  the  ground,  a  feature  not  always  available  with  a 
conventional  RNAV  unit  due  to  VOR  coverage. 

It  was  found  that  the  GPS  systea  was  able  to  maintain  navigation  in  the 
mountainous  area  near  the  Burlington,  VT  Airport  and  in  the  vicinity  of  large 
buildings  at  Logan  International  Airport.  The  systea  was  able  to  maintain 
satellites  in  track  to  elevation  angles  as  low  as  5*,  even  when  on  the 
ground. 

The  systea  was  generally  able  to  maintain  satellites  in  track  during 
turns  with  only  teaporary  loss  of  lock.  The  systea  was  able  to  recover  froa 
the  teaporary  loss  of  satellite  without  a  substantial  effect  on  the  systea 
position  accuracy.  It  was  seen  that  the  horizontal  position  accuracy  was 
aalntained  when  the  systea  was  forced  to  navigate  with  three  satellites  during 
the  teaporary  loss  of  one  satellite  in  a  turn. 


TABLE  4-16 


POSITION  ACCURACY  STATISTICS  FOR  OPERATIONAL  FLIGHT  TESTS 

4-10  FEBRUARY  1983 


Position  Fixes: 

Error 

Mean 

Horizontal 

98 

Vertical 

103 

Tracker  Estiiaates: 

Error 

Mean 

Horizontal 

98 

Vertical 

115 

Standard 


Deviation 

RMS 

95 % 

143 

174 

305  ft 

78 

129 

242  ft 

Standard 


Deviation 

RMS 

952 

161 

189 

333  ft 

93 

148 

277  ft 

229 
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5.0  FUNCTIONAL  REQUIREMENTS  OF  A  GENERAL  AVIATION  RECEIVER 

The  results  of  testing  the  GPS  Tests  and  Evaluation  System  evaluation  of 
the  GPS  T  and  E  equipment  Indicate  that  a  general  aviation  GPS  navigator 
should  satisfy  the  following  functional  requirements: 

5.1  General  Requirements 

5.1.1  Reliability 

The  GPS  navigation  system  encompassing  the  user  equipment,  satellite 
vehicles  and  ground  support  equipment,  must  have  a  combined  reliability  equal 
to  or  surpassing  that  of  alternative  navigation  systems.  The  GPS  user 
equipment  must  therefore  provide  navigation  data  without  operationally 
significant  outages  due  to  fades  or  multipath,  assuming  that  the  NAVSTAR 
constellation  provides  acceptable  coverage. 

5.1.2  Integrity 

The  GPS  navigation  system  must  not  provide  misleading  information  under 
any  conditions  of  operational  significance.  It  is  therefore  necessary  that 
the  GPS  receiver  continually  monitor  its  own  performance  and  indicate  to  the 
pilot  when  the  navigation  information  displayed  is  no  longer  in  compliance 
with  the  accuracy  requirements. 

Further,  provisions  must  be  incorporated  to  allow  the  pilot  to  verify 
that  the  navigation  information  is  accurate  using  either  built-in  test 
equipment,  an  auxiliary  test  system,  or  a  procedural  check. 

5.1.3  Compatibility 

The  GPS  navigation  system  must  provide  a  pilot  Interface  which  is 
compatible  with  existing  air  navigation  systems. 

5.2  Technical  Requirements 

5.2.1  Link  Margin 

The  GPS  receiver  link  margin,  for  all  satellites  above  5°  elevation  angle 
and  with  nominal  effective  radiated  power,  should  be  at  least  3  dB  during 
acquisition  and  at  least  6  dB  once  a  satellite  is  in  track.  These  margins 
assume  the  aircraft  is  in  level  flight. 

5.2.2  Startup 

The  receiver  should  automatically  acquire  and  track  all  visible 
satellites. 


5.2.3  Time-to-First-Fix 


The  receiver  should  automatically  provide  a  first  fix  to  the  navigation 
display  within  6  minutes  of  power  on,  unless  the  set  has  not  been  operated  for 
longer  than  one  month. 

5.2.4  Continuous  navigation  Solutions 

The  GPS  receiver  should  provide  position  estimates  or  related  navigation 
data  derived  from  all  satellites  in  view,  to  the  navigation  display  at  a  1  Hz 
rate. 


5.2.5  Position  Estimation 


The  position  estimation  algorithm  should  substitute  a  synthetic  altitude 
(coast)  if  the  HOOP,  computed  each  fix  Interval,  becomes  unacceptably  high,  or 
if  the  number  of  satellites  drops  temporarily  to  three. 

5.2.6  Performance  Monitor 


The  performance  monitor  should  perform  the  following  tests: 


Test 

Test  satellite  health  data. 

Test  pseudo-range  for 
consistency  with 
previous  measurement 
or  present  prediction. 

Test  position  estimation 
variance. 

Maintain  coast  timer. 


Action 

Delete  unhealthy  satellites. 

Do  not  use  inconsistent  values  in 
the  position  estimator. 


Set  NAV  FLAG  to  notify  pilot  of 
bad  data  if  variance  exceeds 
threshold. 

Set  NAV  FLAG  to  notify  pilot  of  bad 
data  if  coast  time  >  T  where  T  is 
selected  so  that  an  aircraft 
maneuver  will  not  cause  an  error 
greater  than  the  applicable  FAA 
accuracy  requirement. 
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APPENDIX  A 

POSITION  ESTIMATION  ALGORITHM 


The  receiver  software  uses  a  least-squares  fix  which  operates  by  batch 
processing  a  set  of  up  to  ten  GPS  pseudo-range  measurements  that  are 
sequentially  obtained  by  the  receiver.  The  algorithm  produces  a  new  fix  every 
2.2  seconds.  The  algorithm  is  outlined  in  Figs.  A-l  and  A-2. 

The  algorithm  begins  by  referencing  all  10  measurements  to  a  common  time 
by  utilizing  the  pseudo-range  rate  estimates  derived  from  the  receivers  AFC 
loop,  i.e, 

PRk<t r>  -  PMtfc)  +  P*k<tk>*ltr-tk] 

where 


tr  ■  common  reference  time 

tfc  *  is  the  measurement  time  of  the  k-th  satellite  and 

PRR  and  PRR  are  the  pseudo-range  and  pseudo-range  rate  for  the  k-th 
satellite. 

The  common  reference  time  is  selected  to  be  the  mid  point  of  the 
2.2  second  measurement  batch  interval  in  order  to  minimize  prediction  errors. 
Since  the  AFC  loop  has  a  tracking  error  standard  deviation  of  7  Hz  (equivalent 
to  a  range  rate  of  4.1  feet  per  second),  the  measurement  error  contributed  by 
the  translation  process  is  =  4  feet,  one-sigma. 

The  pseudo-range  measurements  are  then  corrected  for  propagation  delay 
effects  using  the  models  shown  in  Table  A-l. 

Referring  to  Figure  A-l,  let  £  ■  (x,y,z)  denote  the  user  position  vector 
in  ECEF  (earth-centered,  earth-fixed)  coordinates  and  let  ££  *  (xi*yi»*i) 
denote  the  position  vector  of  the  1-th  GPS  satellite.  The  pseudo-range 
measurements  are  scalars  given  by: 

*i  ■  l£  -  £il  +  t  (1) 

where  t  is  the  receiver  clock  bias,  normalized  by  the  speed  of  light  to  units 
of  distance.  Equation  (1)  is  expanded  results  in  a  Taylor  series  about  the 

A  A  A 

estimated  user  position,  r,  and  estimated  receiver  clock  time  tr  *  t  +  t. 

The  equation  is  then  linearized  by  neglecting  all  but  first-order  terms: 


9 mi  ^m^ 


2 


TABLE  A-l .  SIGNAL  PROPAGATION  DELAY  COMPENSATION  MODELS 


IONOSPHERIC  [MODEL  NOT  IMPLEMENTED] 


eIONO  -  6  (E)  [d x  +  d2  cos  (2n  ( - )]  (feet) 

24 

9  -  1,  E  -  90°,  0  »  3,  E  -  5° 

E  ■  elevation  angle  to  the  satellite 
dj  and  d2  are  data  base  parameters 
t  Is  local  tine  in  hours 


TROPOSPHERIC 


'TROPO 


K  e“h/hs 


(feet) 


K  -  N0h8 


N0  “  Sea  level  refractivity 
h  “  estimated  altitude 


h8  ■  Exponential  scale  height  (22500  ft  above  mean 
earth  radius) 

h  “  estinated  altitude 

E  ■  satellite  elevation  angle 


\  •*.  i 


yy< 
y  \  y 


:Mx*. 

•>y^ 

•  ■ 


-SV-V* 


where  A_r  -  r  -  r  and  Ab  ■  tr  -  tr  ■  t  -  t.  The  Am^  are  the  residual 
pseudo-range  errors  due  to  the  user  position  estimate  error,  Aj:,  and  receiver 
clock  bias  estimate  error,  Ab.  Rewriting  equation  (1),  we  have: 

■i  ■  [(X“X£>2  +  (y  -  +  (z  -  +  tf  “  * 

thus 

(x  -  x*)  (y  -  y*)  (s  -  Zi> 

Am^  - - Ax  + - Ay  +  — - - -  Az  -  Ab 

\l  -  2±  I  \L  ■  £i  I  IjL  -  £±  I 

Ar 

■  f^il*  ^12*  *4.3*  *  I  1 

Ab 

where  the  hjj  are  the  direction  cosines  from  satellite  1  to  the  user 
position. 

When  N  measurements  are  available  (N  >  4),  Am  in  matrix  form  becomes : 

Ar 

Am  ■  H  •  [  —  ] 

Ar 

-  H  •  AS 

where  H  is  an  N  x  4  matrix  and  Am  is  an  N-dimenslonal  vector.  Solving  for  AS: 
Ar 

AS  -  —  -  (HTH)-1  HT  Am  (3) 

Ab  ” 

In  the  special  case  of  N-4 ,  then  (HTH)“*  HT  ■  IT*. 

The  resulting  values  of  A£  and  Ab  are  used  to  correct  the  previous  user 
position  and  clock  bias  estimates  to  produce  new  estimates  for  the  next 
measurement  cycle. 

The  batch-processing  method  allows  the  data  for  all  satellites  In  view  to 
be  utilized  to  produce  a  position  fix,  and  obviates  the  need  for  a  satellite 
selection  algorithm,  since  the  position  fix  utilizing  all  N  satellites  will 
have  better  GDOP  than  any  subset  of  four  of  the  satellites.  Furthermore,  It 
can  be  shown  that  the  residual  error  is  minimized  in  the  least-squares  sense 
by  equation  (3),  under  the  linearity  assumption  of  equation  (2). 


APPENDIX  B 


POSITION  ESTIMATION  FILTER 


The  position  estimates  developed  by  the  algorithm  described  in  Appendix  A 
are  processed  by  a  filter  (tracker)  before  they  are  released  to  the  navigation 
software. 

An  alpha-beta  fixed  gain  filter  was  selected  following  a  study  of  tracker 

•  • 

alternatives.  The  filter  operates  on  X,  Y,  Z  and  t  to  produce  X.  Y,  Z,  X,  Y, 

•  •  n 

Z,  t  and  t.  The  velocity  gain  is  related  to  the  position  gain  by  B  *  oV( 2-o). 
Figure  B-l  shows  the  general  flow  for  the  tracker. 


FIG.  B-1.  POSITION  ESTIMATION  TRACKER 
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APPENDIX  C 


COORDINATE  CONVERSION  ALGORITHM 

In  order  to  compare  the  actual  aircraft  position,  as  estimated  by  the 
Mode  S  Facility,  with  recorded  GPS  T&E  position  estimates  it  is  necessary  to 
to  convert  the  GPS  estimate  coordinates  to  Mode  S  Facility  local  coordinates. 
The  latter  is  described  in  the  following  paragraphs. 

The  three  coordinate  systems  of  concern  are:  (1)  earth  centered  earth 
fixed  (ECEF) ,  (2)  local,  and  (3)  geodetic.  The  ECEF  system  Is  a  Cartesian 
coordinate  system  fixed  to  the  earth  and  defined  such  that  the  z-axis  points 
toward  the  North  pole,  the  x-axls  points  out  along  the  intersection  of  the 
Greenwich  meridian  and  the  equator,  with  the  y-axis  chosen  to  complete  a 
right-handed  coordinate  system.  A  local  system  is  constructed  so  that  the 
x-axis  points  east,  the  y-axis  points  north  and  the  z-axis  points  "up".  The 
geodetic  coordinates  consist  of  longitude,  latitude,  and  altitude  above  sea 
level. 


Conversions  between  these  coordinate  systems  depend  on  the  model  used  for 
the  surface  of  the  earth.  For  this  effort  the  WGS-72  model  was  selected. 
WGS-72  is  an  ellipsoid  of  revolution  whose  semi-major  (equatorial)  radius  is 
6378135.0  meters,  and  whose  semi -minor  (polar)  radius  is  6356750.5  meters. 

Conversion  from  local  coordinates  to  ECEF  (and  inversely)  requires  a 
rotation  and  a  translation.  (The  rotation,  in  effect,  aligns  the  coordinate 
axas  and  the  translation  displaces  the  origin).  If  P  is  a  point  whose 
coordinates  in  the  local  system  are  expressed  as  the  three-vector  quantity 
v  -  (x,y,s),  then  its  coordinates  in  ECEF,  V  •  (X,Y,Z)  are  given,  in  general, 
by: 


V  -  Rv  +  D 


[1] 


where  R  is  a  3x3  rotation  matrix  and  D  is  the  ECEF  three-vector  coordinate  of 
the  origin  of  the  local  system.  The  Inverse  conversion,  is  given  by 


v  -  Rt  (V  -  D) 


I2J 


where  R?  is  the  transpose  of  R  (which  is  equal  to  its  inverse). 


To  state  the  transformation  between  a  local  system  and  ECEF  coordinates, 
let  6  be  the  longitude,  A  be  the  geodetic  latitude  and  A’  be  the  geocentric 
latitude  of  the  local  coordinate  system.  Then  the  rotation  matrix  is  given 
by: 


R 


-sin  6  -  cos  9  sin  A 

cos  6  -  sin  8  sin  A 

0  cos  A 


cos  6  cos  A 
sin  8  cos  A 
sin  A 


[3] 


C-l 


To  coapute  D,  we  need  to  consider  the  following  geoaetry: 


a 


2  2 
X  Z 

the  nuaerical  eccentricity.  Prom  the  equation  for  an  elipse,  — 2  +  — 2  -  1 

a  b 


it  can  be  shown  that: 

D-  ■  (d(l-e2)  +  h]  sin  X 

[4] 

Dx  *  (d  +  h)  cos  X  cos  6 
Dy  ■  (d  +  h)  cos  X  sin  0 

where: 


a 

d  -  — - 

/l-e2  8in2X 

and  h  is  the  height  (altitude)  of  the  local  systea.  Thus  equations  [3]  and 
[4]  convert  local  coordinates  to  ECEF  coordinates. 

Using  R  and  D,  conversion  froa  local  coordinates  to  ECEF  coordinates  and 
back  can  be  accomplished.  If,  for  exaaple,  the  coordinates  of  a  satellite  are 
given  in  a  horizon  systea  (azimuth,  elevation,  and  slant  range)  these 
coordinates  can  be  converted  to  local  x,y ,z  and  then  to  ECEF.  A  difficulty 
arises,  however,  when  converting  the  Mode  S  Facility  reported  position  of  a 
target  to  its  own  local  coordinates.  The  difficulty  arises  because  target 
reports  give  altitude  rather  than  elevation  angle.  Because  of  this,  the  local 
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